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SECTION 1.0 


SUMMARY 


In an examination of component combustor noise data it was found that a 
sound power level correlation could be arrived at based on the steady state 
heat release rate and inlet temperature. This correlation is presently valid 
only for component rig and the configurations tested.. Before this correlation 
can be applied to engine data, the effect of blade row attenuation and exhaust 
nozzle termination must be defined. 

The sound pressure levels and spectral shapes of too combustors were 
recorded. The CF6-50 was the baseline combustor. This combustor design is 
currently being used in the CF6-50 engine. The Double Annular Combustor 
configuration D13 was one of the final clean combustor design configurations 
tested under Phase II. The two combustors were tested at pressures up to 
9.5 atmospheres. Phase I testing was conducted at pressures of 4.8 atmo- 
spheres or less. Only configuration D13 was tested at both the higher Phase 
II pressures and the lower Phase I levels. A combustor noise prediction 
correlation was arrived at based on this configuration and was found to be in 
good agreement with the other configurations. 

The analysis of noise attenuation due to turbine blade rows was extended 
to include the supersonic relative flow condition. The predicted attenuation 
based on this analysis was calculated to be greater than 34dB for the CF6-50. 
However, this level of attenuation has not been verified by test data. Work 
in this area has been initiated under NASA contract number NAS 3- 19 4 35 . 



SECTION 2.0 


INTRODUCTION 


Low frequency noise radiating from the core exhaust of gas turbine 
engines contributes to an engine’s overall perceived noise level, particularly 
at approach. The lower frequency core noise is related to the engine’s 
combustor. In order to predict and possibly reduce combustor noise at the 
source It was necessary to determine the controlling parameters. This was 
accomplished utilizing component combustor test data. A relationship was 
identified which correlated the overall combustor noise power levels to 
easily determinable aerodynamic and thermodynamic parameters. This empirical 
relationship may make it possible to extrapolate the noise level for an 
annular combustor in a component test rig to actual engine conditions. It 
remains to be seen whether this equation is valid for a combustor installed 
in an engine. 

In calculating the overall noise levels it was assumed that the measured 
pressure fluctuations in the combustor were dominated by acoustic waves 
rather than convected disturbances. This assumption was believed to be 
correct, however it remains to he verified. 

The information contained here in Phase II is a continuation of the 
work presented under Phase I of this program. The purpose of this phase was 
to 1) measure noise from two combustors, 2) correlate these and the data 
from Phase I with their respective aerodynamic and thermodynamic parameters 
and 3) analyze the effect of turbine blade row attenuation. Two combustors 
were chosen for testing. The CF6-50 was selected as a current design com- 
bustor. This permitted data taken from this combustor to be evaluated 
against farfield noise levels from the CF6-50 engine. The Double Annular 
configuration D13 was selected as an advanced design clean combustor. The 
data from this configuration was compared to the current design in order to 
determine the effect of the design change on the noise. 

This report supplies supplemental information to the Experimental Clean 
Combustor Program Phase II Final Report (Reference 1). The test point data 
and the corresponding reading numbers in this Noise Addendum and in Reference 
1 are the same. 
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SECTION 3.0 


INSTRUMENTATION 


Both the aerodynamic and acoustic instrumentation were located at two 
axial planes in the component combustor test rig (see Figure 1). The upstream 
location was designated engine station 3.0 which is typically the compressor 
discharge plane in an engine (see Figure 2) . The downstream instrumentation 
was located directly behind the combustor and was designated engine station 
3.9. Reference to a station 4.0 refers to a plane immediately downstream of 
the combustor where the turbine first stage nozzle would be in an engine. 

Acoustic data were taken utilizing waveguide probes and kulite sensors. 
The upstream acoustic waveguide probe was 50.8 cm long and was located at 
45° F.L.A. (Forward Looking Aft) from the top centerline. The downstream 
acoustic probe sensing holes were located at the combustion exit. For the 
Double Annular Combustor test a wall static probe 35.56 cm long was located 
at 84° F.L.A. This probe was flush with the inside of the deflector plate 
and was axially located in the same plane as the body of the downstream 
waveguide probe. 

The aerodynamic data were taken using standard instrumentation available 
in the CF6-50 component combustor rig. Information regarding details of the 
aerodynamic instrumentation of the test facility can be found in Reference 1. 


2 , 1 ACOUSTIC WAVEGUIDE PROBE DESIGN 

Due to the high temperature environment of a combustor the acoustic 
data were acquired through the use of waveguide probes. This permitted the 
pressure transducers to be located in a relatively lower temperature area. 

A waveguide probe consisted of a hollow tube with three parts: the sensing 

end, the transducer, and the infinite or pc termination. The probe functioned 
by sensing pressures at the sensing end and permitted the pressure to propagate 
through the .waveguide tube, past the transducer, to be dissipated without 
reflection in the infinite or pc termination. The sensing end may be either 
a wall static pressure tap installed with the end termination flush with a 
test duct wall, or may extend into the flow in the duct. The probe design 
whose sensing end extends into the flow from a side wall was bent 90 degrees 
into the direction of the flow and a bullet nose tip was inserted in it’s 
end. Pressure entered the waveguide through 36 radial holes immediately 
behind the tip. The sum of the areas of the sensing holes corresponded to 
the cross sectional area of the waveguide, thereby minimizing any pressure 
loss due to area changes. From the sensing end the pressure waves passed 
through the waveguide tube past the transducer. The Kulite transducer liras 
mounted in a compression fitting tee so that the diaphram was essentially 
flush with the inside wall of the waveguide. The Kulite that was used for 
recording data was model XGE-1S- 3 7 5V-200D . Once past the transducer the 
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pressure waves were dissipated by air viscous loss in a continuation of the 
waveguide, which was over 12 meters long. This long length of waveguide was 
called the infinite or pc termination because waves could not be reflected 
back from it. Reflected waves had to be eliminated because they would 
generate standing waves through the waveguide which would generate mislead- 
ing levels at the transducer. 

The waveguides consisted of 6.35 mm diameter tubes with 0.71 mm wall 
thickness. The upstream probe was made of stainless steel and was capable 
of operating at temperatures up to 1000° K, pressures up to 200 atmospheres 
and Mach numbers over 0.5. The downstream probe was sheathed in a water 
jacket which incorporated tip impingement cooling. Details of the downstream 
probe are show in Figure 3. This probe was capable of operating at tempera- 
tures over 1600° K, pressures up to 200 atmospheres and Mach numbers of 0.4. 
In order to meet that condition the cooling water pressure differential 
through the probe was kept at 6,8 atmospheres. At that pressure the water 
flow rate was four liters per minute. 

3.2 ACOUSTIC WAVEGUIDE PROBE CALIBRATION 

The waveguide probe with the water cooled casing was tested in a wave- 
guide probe calibration facility (Figure 4). The calibration consisted of 
measuring the sound pressure level difference between the probes sensing 
holes and the transducer end. A small (1.27 mm dia.) reference Kulite was 
temporarily attached to the probe tip next to the static sensing holes. The 
location was such that the phase associated with the distance between the 
sensing holes and the Kulite was less than nine percent of the wavelength at 
5000 Ha. The percent of phase difference becomes smaller at lower fre- 
quencies. The reference Kulite had a flat response from 20 to over 10,000 
Ha. A B&K type 4134 microphone was connected to the probe's transducer end.' 
The probe tip and reference Kulite were then placed in the plane wave tube 
which was excited at one end by a loudspeaker driver. The speaker was used 
to generate pure tone levels from 50 to over 5000 Hz. In order to keep a 
constant sound pressure level at the probe sensing holes, the reference 
Kulite was connected to a feed back compressor circuit. The signal from 
both the Kulite and the microphone were filtered with a 20 Hz constant 
bandwidth filter centered around the speaker driver frequency. The measured 
difference in level between the two transducers was the probe loss, see 
Figure 5. The calculated loss also shown on the plot was based on the 
Iberall method (Ref. 2). It was observed that the measured response and 
calculated response in general were within one dB. It was therefore con- 
cluded that the effect of the small transducer near the sensing holes was 
negligible. The tabulation of probe attenuations used in correcting the 
one-third octave band data is listed in Table I. 

3.3 EFFECT OF FLOW ON ACOUSTIC PROBES 

Since the probes were to be immersed in the combustor flow, a test was 
conducted in order to evaluate the impact of any probe-generated noise on 
the measuring transducer. The waveguide probes were mounted in a 12.7 by 
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Figure 3. Water Cooled Acoustic Waveguide Probe. 






















Table I. Tabulation of Waveguide Probe Attenuation. 


Frequency 

Upstream 

Waveguide 

Probe 

Downstream 
Water-Cooled 
Waveguide Probe 

— ■ — — - 

Downstream 
Wall Static 
Waveguide 

31.5 

0,30 

2.0 

0.21 

40 

0.34 

2.1 

0.24 

50 

0.38 

2.1 

0.27 

53 

0.43 

2.1 

0.30 

SO 

0.4S 

2.2 

0.34 

100 

0.54 

2.3 

0.38 

125 

0.60 

2.3 

0.42 

160 

0.68 

2.4 

0.47 

200 

0.76 

2.5 

0,53 

250 

0.85 

2.6 

0..60 

315 

0.96 

2.7 

0.67 

400 

1.07 

2.9 

0.75 

500 

1.20 

3.0 

0.84 

630 

1.35 

3.2 

0.94 

800 

1.51 

3.4 

1.06 

1000 

1.70 

3.6 

1.19 

1250 

1.91 

3.9 

1.33 

1600 

2,14 

4.3 

1.50 

2000 

2.40 

4.8 

1.68 

2500 

2.69 

5.3 

1.88 

3150 

3.02 

6.0 

2 . 11 

4000 

3.39 

6.7 

2.37 

5000 

3.80 

7.3 

2.66 

6300 ' 

4.27 

8.3 

2.99 

8000 

4.79 

9. 2 

3.35 

10000 

5.37 

10.3 

3.79 

Frequency in Hz, Attenuation in dB 




6.1 cm duct, which simulated the combustor annulus passage height. The duct 
had a muffled upstream plenum and three flow-straightening screens, see 
Figure 6. From the plenum to the duct the flow was accelerated, in order to 
generate a thin boundary layer. The measured centerline turbulence levels 
were low, on the order of 1.1%. The noise level of this facility was also 
low. The overall sound pressure level measured at the wall without the 
probe installed was 120 dB at a Mach number of 0.3. 

The two probes that were used in the combustor noise test were each in 
turn mounted in the duct. The sound pressure level spectra for each were 
measured over the range of Mach numbers from 0.1 to 0.3, Figure 7. At the 
0.3 Mach number the 1.27 cm diameter downstream probe had an overall sound 
pressure level of 129.3 dB while the 0.635 cm diameter upstream probe had 
132 dB 0ASPL. The difference in level between the two probes was probably 
due to the difference in the shape of the probe tips. The flow-generated 
noise over the probes was on the order of 130 dB, ten dB more than the 
facility noise. Also, the overall pressure levels measured in the combustor 
test were all over 145 dB. The measured combustor overall pressure levels 
were therefore 15 dB or more over probe generated noise levels. It appeared 
that probe flow noise would not influence the combustor noise measurements. 

The effect of water flowing through the passages in the downstream 
probe wauer jacket was also tested. Figure 8(a) shows that the effect of 
the water flowing through the probe was negligible when compared to the 
water off. The peak level was also very low, less than 80 dB. This test 
was performed while the probe was still installed in the flow noise duct. A 
peak in the spectra was observed at 315 Hz, which also occurred in the flow 
noise spectra. The peak was not, however observed with the wall microphone. 
Furthermore a test (Volume II of Reference 4) was conducted with nothing in 
the duct. An average of eight microphones located in the room around the 
duct and a flow Mach number of 0.3 produced similar peaks. Figure 8(b), at 
315 and 2500 Hz which were comparable to Figure 7. The peaks therefore 
appeared to be standing waves related to the duct. 


3.4 THE ACOUSTIC MEASUREMENT SYSTEM 

The electrical system (Figure 9) consisted of Kulites connected to low 
noise Princeton Applied Research amplifiers which fed into a Lockheed wide- 
band FM tape recorder. The Kulites were Model XCEL-31-13-10QD which, have a 
flat frequency response from 20 to over 10,000 Hz. The wall static probe 
that was used in the Double Annular combustor test as well as the waveguide 
probes were all operated in this manner. 

A backflow purge was used on all of the probes. The purge consisted of 
connecting the end of the pc termination to a higher static pressure upstream 
of the combustor. The flow rate was restricted with a 3.18 mm diameter 
tube, 30.5 cm long. A back pressure of less than 0. 9 atmospheres produced 
less than 24.8 liters per minute air flow for two probes. The peak 1/3- 
octave levels at this condition was less than 110 dB. This level was con- 
sidered negligible for this test. The actual pressure across the system 




Figure 6, Acoustic Duct Facility for Probe Flow Noise Mesurentent. 
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during the combustor tests varied with air flow rate. However, it did not 
exceed the 0,9 atmosphere level. Therefore the purge noise levels would be 
even lower. This system was the same as the one used under Phase I (Reference 
3) whose spectral shape and levels can be seen in Figure 10. 


16 
















SECTION 4.0 


COMB'D STOR DESIGNS 


4.1 THE CF6-50 ENGINE 

The CF6-50 engine is one of too versions of the CF6 high bypass ratio 
turbofan engine. The other model, the CF6-6D, is the power plant for the 
McDonnell-Douglas DC-10 Series 10 Tri-Jet aircraft. The CF6-50 is the higher 
power version of the two engines and is being utilized in the McDonnell- 
Douglas DC-10 Series 30 Tri-Jet long range intercontinental aircraft, the 
Boeing 747, and the Airbus Industrie A300B aircraft. Some of the overall 
specifications of the CF6-50A engine are: 


Weight 

Length (cold) 

Max. Dia. (cold) 
Fan/Comp. Stages 
HPT/LPT Stages 
Thrust /Weight 
Pressure Ratio 
Airflow 

Max . SLS Thrust 
SFC 

Cruise Mach Number /Alt 

Thrust 

SFC 

Bypass Weight flow Ratio 


3780 kg 
482 cm 
272 cm 
l-(3)/14 
2/4 
5.95 
30:1 
650 kg/s 
218 kN 
0.389 

0. 85/10.5 km 
48 kN 
0.654 

(Fan/Core) 4.1 


The smoke emission levels of this engine are low and virtually invisible 
at all operating conditions. The favorable fuel economy and high thrust— to— 
weight ratio made this engine an attractive selection as a baseline vehicle. 


4,2 CF6-50 COMBUSTOR 


f The CF6-50 combustor is a high performance combustor incorporating a 

■' low pressure loss step diffuser, carbureting swirl— cup-dome design and a short 
burning length. The step diffuser has a pressure loss of about one percent 
of the total pressure and does not vary significantly over the engine cycle. 
These features have made it possible for this combustor to have low smoke 
emissions and high combustion efficiency at all operating conditions, a low 
pressure loss, and low exit temperature pattern factors. The design of the 
step diffuser provides a uniform, steady airflow distribution into the com- 
bustor. Its shorter burning length which requires less liner cooling air 
improves the exit temperature pattern factor and profile. 


IS 



The flow at the compressor discharge plane has a Mach number of 0.27. 

This high velocity flow is diffused in a relatively long, area-rule predif- 
fuser with an area ratio of 2.0. The top half cross section. Figure 11(a), 
illustrates the combustion system layout. The production CF6-50 combustor 
incorporates 30 vortex- inducing axial swirl er cups, one for each of its 
corresponding fuel nosales. 

The air that flowed through the prediffuser passage was W 3 . The air 
which passed through the combustor exit was W c . The difference between the 
two was the simulated turbine bleed air. Past the prediffuser passage the 
air flow was split into three streams, outer, center and inner. The com- 
bustion air flow W c entered the combustion sone through either swirl cups, 
dilution holes, or film cooling slots. The outer and inner flow streams 
were smoothly accelerated around the combustor cowling and entered the com- 
bustion zone through the outer and inner liner passages which consisted of 
holes and slots. Thirty percent of W c was liner cooling flow which centered 
the combustion reaction zone through film slots. The center stream flowed 
into the combustor primary zone through the dome cowling to the swirl cups. 

The cowling opening was oversized and provided free-stream diffusion of the 
dome flow. This resulted in higher pressure recovery in the center stream. 

The higher pressure recovery generated higher pressure drops across the dome 
and therefore higher velocities through the swirl cups and other dome flow 
openings. 

The combustion process of chemically reacting a liquid fuel with air, 
for heat release, is by nature nonhomogeneous. The local fuel-air ratio will 
be a function of the location in the combustor. When the fuel enters the 
primary zone close to the swirl cups, the fuel-air ratio is close to stoichio- 
metric and decreases from that value as you travel axially down the combustor. 
The dillution holes and film cooling slots dilute the primary zone mixture, 
and to an extent react with it. Figure 11(b) shows the details of the air 
flow splits. It may be useful to determine this spatial function of the 
fuel-air ratio in order to assess its effect on the noise. However, in order 
to simplify this study of combustor noise, the fuel-air ratio that passed 
through the combustor exit was used. 

Relevant design parameters for the CF6-50 combustor are shown in Table 
II. The space rate listed in the table is the heat released in the combustor. 
It is defined as the actual heat release, in joules per hour, divided by the 
product of the total enclosed combustor volume down to Station 4 (in cubic 
meters) and the burner total pressure in atmospheres. 


Space rate = 


heat release 
volume * pressure 





Table II, CF6-5Q Combustor Design Parameters, 


Combustor Airflow (W c ) 
Compressor Exit Mach Number 
Overall System Length 
Burning Length (L B ) 

Dome Height (Hj,) 

l b /h d 

Reference Velocity 
Reference Area 
Space Rate 
AP t /Pt 3 

Number of Fuel Nozzles 
Fuel Nozzle Spacing (B) 

l b /b 

b/h d 


103.42 kg/s 
0.27 

75.95 cm 

34.8 cm 
11.43 cm 

3.0 

25.9 m/s 
3729 cm 2 

2.2 x 10 11 J/hr-m 3 
4.3% (Total) 

3° 

6.91 cm 

5.0 


-atm 


0.60 



4.3 DOUBLE ANNULAR COMBUSTOR 

The general arrangement of the Double Annular Combustor design approach 
is shown in Figures 12 and 13 and its geometric parameters are tabulated in 
Table III. The dome assembly consisted of two annular spectacle plates 
separated by a small centerbody. Assembled in the spectacle plates was an 
array of 60 air swirlers (30 in each annulus) . The air swirler components 
consisted of a primary air swirler /venturi casting, a counter-rotating 
secondary air swirler, a flameshield and a retainer ring. The air swirler 
assemblies were attached to the dome spectacle plates with a radial slip 
joint arrangement to accommodate mechanical stack-up and thermal growth 
between the fuel Injectors and the combustor assembly. The flameshields were 
impingement cooled. The centerbody and dome panels were film cooled using 
wiggle strip construction. 

The cooling liners x*ere modified production CF6-50 liners. The cooling 
flow metering hole areas were reduced and the dilution holes were modified, 
Figure 14, 
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SECTION 5.0 


EXPERIMENTAL RESULTS 


The results presented here consist of tables of the aerodynamic data 
and one-third octave band sound pressure levels. Plots of the sound pressure 
spectra are also presented. 

5.1 AERODYNAMIC DATA 

The aerodynamic data consist of the fuel-air ratio, total and static 
pressures and total temperatures at both upstream and downstream locations 
in the combustor, the air and fuel flow rates and the fuel temperature for 
each reading number. A total of eight data points were recorded for the 
CF6-50 combustor, four fuel-air ratios at both the approach and takeoff inlet 
conditions. The dry loss, zero fuel air ratio test points could not be taken 
due to a cokeing problem with the fuel nozzles. A low fuel air ratio was 
substituted. The aerodynamic data are listed in Table IV. The combustor 
pressure levels were higher than those tested under Phase I. Thus a one to 
one comparison of this combustor to any of the Phase X configurations is not 
recommended. 

The Double Annular aerodynamic data are listed in Table V. A total of 
eighteen data points were recorded encompassing approach and takeoff inlet 
conditions comparable to those taken for both the CF6-50 and the Phase I 
combustors. The 3.4 and 4.7 atm pressures correspond to Phase I approach 
and takeoff while the 6.8 and 9.5 atm pressures correspond to the CF6-50 
Phase II approach and takeoff pressures respectively. The aerodynamic data 
listed in Table V are presented in terms of successive reading numbers with 
the exception of readings 748 and 699. Those two readings were included in 
their group of similar inlet conditions. This order is maintained consistent 
through the reporting of the acoustic data-tabulations and their plots. 


5.2 ACOUSTIC DATA 


The one-third octave band sound pressure levels from 31.5 to 10,000 Hz 
and the calculated overall sound pressure level based on that spectrum are 
listed in Tables VI and VII for the CF6-50 and Double Annular combustors 
respectively. All of the levels have been corrected for probe losses. The 
notation U, indicates the upstream probe, D the downstream probe and the sub- 
script W indicaf s the wall static probe. 

During the Double Annular combustor test it was not possible to' get data 
from the upstream probe for all of the readings. However, data taken with a 
wall static probe at the downstream location have been substituted. Reliable 
data from the upstream probe were obtained for readings 748, 723, 726, 727 
and 728. The CF6-50 combustor noise spectra are plotted on Figures 15 to 22. 
The Double Annular combustor were plotted on Figures 23 to 45. 
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Table IV. CV6-50 Aerodynamic Data 




(atm. 

) 

(kgm/sec) 

— 11 

<° 10 

KDG. 

f /a 

P T3 

P S3 

P T3.9 

P S3.9 


W f 

t T3 


■•'fuel 

317 

.0141 

6.87 

6.41 

6.45 

6.29 

29.9 

1521 

661 

1174 

295 

318 

.0169 

6.80 

6.40 

6.43 

6.27 

30.4 

1852 

659 

1266 

295 

319 

.0197 

6.87 

6.48 

6.51 

6.33 

30.6 

2171 

659 

1358 

295 

320 

.0240 

7,28 

6, 86 

6.88 

6.71 

30.2 

2607 

659 

1489 

294 

321 

.0229 

9.59 

9.11 

9.16 

9.00 

33.1 

. 

. 

2733 

811 

1583 

294 

322 

.0199 

9.53 

9,10 

9 .15 

9.00 

33.0 

2361 

818 

1496 

295 

323 

.0169 

9.59 

9.11 

9.15 

9.01 

32.9 

1999 

819 

1402 

295 

324 

.0138 

9.53 

9.09 

9.14 

9.00 

33.2 

1651 

819 

1304 

296 



















Table V, Double Annular, D-13, Aerodynamic Data 




(atm.) 

(kgm/sec) 

(legm/hr) 

no 1 

RDG. 

f /a 

P T3 

P S3 

P T3.9 

P S3,9 

w c 

Wf 

T T3 

T T3.9 

Tr 

705 

.0118 

6.84 

6.50 

6.48 

6.37 

28.0 

1186 

627 

1064 

304 

706 

.0098 

6.82 

6.48 

6.46 

6.36 

28.0 

997 

627 

995 

304 

707 

.0078 

6.83 

6.48 

6.48 

6.38 

28.1 

791 

627 

924 

304 

708 

.0059 

6.84 

6.50 

6.48 

6.39 

27.9 

590 

627 

852 

304 

709 

.0136 

4.84 

4.60 

4.59 

4.52 

17.0 

833 

816 

1291 

302 

710 

.0170 

4.77 

4.54 

4.52 

4.45 

16.6 

1017 

815 

1401 

301 

711 

.0200 

4.77 

4.55 

4.51 

4.43 

16.6 

1195 

821 

1502 

301 

748 

.0234 

. 4.72 

4.51 

4.47 

4.38 

16.5 

1387 

820 

1604 

302 

712 

.0140 

3.45 

3.29 

3.28 

3.22 

13.9 

701 

630 

1139 

302 

713 

.0122 

3.43 

3.26 

3.26 

3.21 

13.7 

602 

629 

1077 

301 

714 

.0101 

3.44 

3.28 

3.27 

3.22 

13.8 

499 

628 

1004 

301 

715 

.0081 

3.43 

3.27 

3.25 

3.20 

13.8 

401 

632 

9 37 

301 

716 

.0061 

3.42 

3.26 

3,25 

3.20 

13.8 

304 

631 

865 

301 

699 

.0229 

9.55 

9.0* 

9.05 

8.88 

33.4 

2752 

819 

1587 

30 3 

723 

.0230 

9.53 

9.07 

9.01 

8.84 

33.3 

2756 

816 

1589 

304 

726 

.0199 

9.55 

9.14 

9.04 

8.88 

33.1 

2373 

820 

1499 

304 

727 

.0170 

9.54 

9.07 

9.04 

8.89 

33.1 

2021 

820 

1405 

304 

728 

.0140 

9.53 

9.07 

9.02 

8.88 

33.1 

1667 

820 

1306 

303 | 
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Table VI. One-Third Octave Band Sound Pressure Levels, CF6-50 Combustor. 



Rdg. 

317 

318 

319 

320 

321 


322 


323 


324 

Probe 

U 

D 

0 

D 

U 

1) 

0 

D 

11 

D 

U 

D 

U 

D 

U 

B 

Frequency 

















31.5 

131.3 

131.5 

131.0 

133.0 

131.0 

133.3 

13D.3 

133.3 

132.5 


132.9 

135.3 

132.0 


132.3 

135,5 

40 


131.2 


131.7 

130.4 

132.0 

129.9 

132.2 

132.9 


132.1 

135.5 

132.4 


132.1 

135.7 

50 


132,1 

130.8 

133.1 

130.8 

132.1 

129.8 

132,9 

133.5 


133.0 

137.6 

133.3 

137.9 

133.0 

137.4 

63 

130.8 

133.1 

130.8 

134.4 

131.0 

134.1 

131.3 

135.4 

134,0 

140.4 

133.5 

138.9 

133.3 

138.9 

133.5 

138.1 

80 

nam 

134.6 

131.1 

134.6 

131.1 

134.6 

130.6 

135.3 

132.4 

139.8 

132.4 

13B.1 

131.9 

138.6 

131.9 

137.3 

100 

132.9 

135.2 

SI 

135.7 

133.4 

135.9 

132.4 

136.7 

133.9 

139.4 

133.1 

137.7 

133.4 

137.9 

133.4 

137.2 

1Z5 

135.8 

136.0 

PK if !■ 

135.3 

135,8 

134.5 


135.3 

137.5 

137.0 

137.8 

136.3 

137.8 

136.0 

137.8 

136.8 

160 

135.0 

138.3 

Kit il ifl 

138.0 

135.2 

13B.3 


138.5 

134.5 

138.5 

134.7 

137.5 

134.7 

138 0 

135.0 

138.5 

200 

137.0 

138.3 

136.7 

137.3 

136.7 

136.8 

136.0 

137.0 

138.0 

136. B 

138.2 

136.0 

138.2 

135,5 

138.0 

136.0 

250 

134.6 

139.1 

134.1 

138.1 

133.9 

137.9 

134.6 

137.9 

134.6 

138.1 

134.6 

137.6 

134.6 

137.6 

134.1 

137.1 

315 

140.5 

140.8 

139.8 

139.0 

140,0 

138.5 

140.3 

140.0 

139.5 

140.0 

138.3 

138.5 

138.3 

138.3 

138.0 

138.3 

400 

146.1 

146.3 

146.6 

146.5 

146.1 

145.5 

149.1 

149.3 

146.1 

143.8 

145.3 

143.0 

145.3 

142.8 

145.1 

143.3 

500 

152.0 

151.7 

152.5 

153.2 

152.5 

153.0 

157.3 

158.0 

151.3 

149.5 

152.3 

148.7 

152.0 

148.7 

152.3 

148.2 

630 

151.0 

149.7 

151.2 

150.2 

150.5 

150.4 

148.7 

148.7 

148.7 

148.2 

148.0 

146.4 

147.7 

146.7 

147.2 

146.4 

800 

151.8 

149.4 

152.5 

149.6 

151.8 

148.1 

149.5 

147. 1 

153.3 

149.9 

153.0 

148.1 

153.0 

148.4 

152.5 

148,4 

1000 

148.0 

147.6 

148.0 

147.3 

148.5 

147.1 

140.2 

147.1 

151.7 

149.1 

150.7 

147.3 

151.0 

147.3 

150.5 

147.3 

1250 

145.0 

149.4 

144.6 

149.1 

145.0 

149.1 

145.5 

148.4 

150.3 

150,1 

149,8 

148.9 

149.5 

149.1 

149.5 

149.1 

1600 

142.8 

147.1 

142.8 

147.6 

142.8 

147.4 

144.1 

145.9 

144.1 

147.4 

143.1 

146 . 4 

143.1 

146.4 

142.8 

146.1 

2000 

140,4 

146.0 

141,1 

145.8 

140.4 

144.5 

141.6 

143.3 

143.4 

145.8 

142.1 

145.0 

142.4 

145.5 

141.9 

145.5 

2500 

137.0 

145.7 

137.8 

145.7 

137.0 

144.7 

138.8 

143.2 

140.5 

146.9 

139 .0 

145.9 

138.8 

146.2 

138.8 

146.4 

3150 

136.2 

144.0 

136.4 

144.0 

136.2 

143.2 

135.9 

141.7 

138.9 

146.0 

137.9 

145.5 

137.4 

145.7 

138.2 

146.2 

4000 

134.1 

141.7 

134.1 

142.0 

133.9 

141.2 

132.9 

140. C 

137.6 

143.5 

138.1 

142.7 

138.1 

142.7 

138.9 

144.2 

5000 

129.6 

140.1 

129.6 

140,3 

129,8 

139. B 

129.3 

139.1 

1 J 3. 3 

143.3 

134.0 

141.3 

134.0 

141.8 

134.8 

142.1 

6300 

121,8 

133.2 

122.8 

133.7 

122.8 

133.5 

122.8 

132.2 

126,3 

ijS „ 7 

126.8 

133.5 

126.1 

134.0 

126. B 

133.2 

8000 

114.9 

126.2 

115.1 

126.7 

115,1 

125.4 

116.6 

125.2 

119.6 

129.0 

118.4 

124.9 

117.9 

125.7 

118.1 

125.4 

10000 

108.1 

121.2 

109.9 

120.7 

108.9 

119.2 

112.1 

119.4 

113.6 

123.9 

113.4 

119.2 

112.6 

119.2 

113.1 

120.7 

OASPL 

15 B. 1 

161.4 

158.5 

159.1 

158.2 

158.7 

159.9 

160.5 

159.2 

158.8 

158.9 

157.6 

158.8 

157.8 

158.6 

157.8 

























RDG. 

705 

706 

i 

70j 

f 

Probe ' 

D 

1 

Dw 

D 

Dw 

D 

Dw 

Frequency 







31.5 

137.3 

135.5 

137.6 

135.2 

137.6 


40 

139.3 

134.6 

139.8 

134.6 

140.0 


50 

140.7 

140.3 i 

141.2 

140.1 : 

141.7 

140.1 

63 

142.2 

137.8 

141.9 

137.5 : 

141.7 

137.3 

80 

140.9 

138.4 

140.6 

137.9 

141.1 

137.4 

100 

139.0 

135.3 

139.5 

135.5 ! 

139.0 

135.0 

125 

139.1 

135.4 

139.6 

135.4 : 

139.1 

134.9 

160 

138.1 

135.9 

137.9 

135.4 

137.1 

134.6 

200 

138.8 

137.3 

138.3 

136.8 

137.3 

135.5 

250 

138,9 

137.6 

137.9 

137.4 

136.2 

135.9 

315 

140.5 

138,9 

140.3 

138.4 

139.0 

136.7 

400 ; 

144.1 

141.4 

145.4 

141.4 

142.4 

139.1 

500 

144.6 

143.6 

147.1 

145.9 

144.1 

143.9 

630 

148.2 

144.1 

151.7 

146.6 

154.2 

147.1 

800 

146.2 

145.6 

147.5 

145.8 

150.5 

147.8 

1000 

148.4 

145.4 

148.6 

I4e>. 1 

148.1 

145.1 

1250 

147.2 

143.3 

149.0 

145.6 

148.0 

145.1 

1600 

145.5 

143.6 

146,2 

145.4 

147.2 

146.1 

2000 

145.8 

141.9 

147.6 

143.1 

149.1 

143.9 

2500 

146.3 

139.5 

148.8 

140,0 

149.8 

141.3 

3150 

145.5 

137.8 

147.0 

140.6 

149.3 

143.3 

4000 

141.3 

134.6 

144.0 

137.1 

146.5 

138.9 

5000 

138.2 

130.7 

139.4 

131.5 

141.7 

131.7 

6300 

135 .6 

129.7 

135.6 

130.7 

136.6 

129.5 

8000 

137.1 

130.0 

137.1 

130.5 

138.4 

129.7 

10000 

140.8 

129.7 

141.8 

130.5 j 

141.6 

: 129.7 


159.2 


160.2 


OASPL 


157.6 


154.4 


155.5 


155.7 


jTiTnr 


E 


708 

709 

710 

D 

Dw 

D 

Dw 

D 

Dw 

138.1 

135.5 

136.1 

132.5 

135.8 

133.0 

139.5 

135.3 

139.0 

133.6 

138.3 

136.8 

140.7 

139.8 

139.7 

136.1 

139.5 

137.6 

140.4 

137.3 

139.4 

133.0 

139.7 

137.8 

139.9 

137.4 

139.4 

132.2 

138.4 

139.7 

138.8 

135.3 

139.0 

132.3 

136.5 

137.0 

137.1 

135.6 

136.8 

132.9 

136.1 

136.4 

136.6 

134.6 

134.1 

133.9 

134.4 

136.6 

136.6 

137.0 

133.6 

134.8 

134.1 

138.0 

135.4 

136.9 

133.9 

135.9 

135.2 

138.9 

137.5 

137.2 

136.8 

135.7 

137.0 

139.7 

140.6 

138.9 

139.6 

136.6 

141.1 

140.9 

140.1 

141.4 

137.6 

138.9 

140.6 

142.9 

146.2 

142.4 

143.0 

139.9 

142.5 

144.1 

146.7 

145.8 

142.5 

142.3 

143.2 

145.1 

147.1 

145.1 

143.9 

140.6 

144.1 

143.6 

144.2 

143.1 

143.2 

141.1 

145 .2 

143.8 

144.0 

144.9 

142.7 

141.6 

143.0 

143.6 

146.3 

141.9 

142.3 

141.6 

142.3 

142.9 

147.0 

140.0 

143.3 

139.8 

143.5 

141.3 

147.8 

143.6 

145.0 

140.1 

145.5 

140.6 

146.5 

141.4 

146.5 

141.9 

146.3 

141,9 

141.4 

134.0 

147.2 

139,2 

147.7 

139.7 

135.1 

130.0 

142.9 

130.2 

143.4 

134.7 

135.9 

129.7 

139.1 

126.2 

139.1 

131.2 

137.6 

127.7 

140.3 

124.0 

139.8 

129.2 

157.1 

154.5 

155.8 

152.2 

156.2 

154.8 
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Table VII. One-Third Octave Band Sound Pressure Levels, 



Frequency 


134. B 

137.5 

138.7 

139.4 

138.1 

137.3 

135.8 

135.1 

134.3 
136.7 

138.3 

141.9 

141.6 

144.0 

145.0 

145.6 

145.2 

143.5 

142.6 

142.5 

144.5 

145.3 

146.4 

141.4 

137.1 


132.7 

136.1 

136.8 

136.8 

137.9 

136.5 

135.6 

135.9 

137.0 

138.6 

139.2 

141.1 

142.4 

143.6 

145.3 

144.9 

143.6 

142.9 

142.9 

141.3 

140.1 

140.4 

139.0 

132.7 

130.5 


129.0 

129.1 

129.4 

129.6 

129.0 

130.4 

133.3 

131.3 

135.5 

134.3 

140.2 

146.2 

147.2 

149.0 

151.5 

151.6 

148.1 

143.4 

143.6 

143.3 

138.7 

134.3 

134.3 

129.7 

131.1 


129.1 

129.5 

130.0 

132.4 

131.9 

131.5 

131.6 

135.6 

135.1 

138.2 

140.3 

142.4 

143.4 

146.0 

146.7 

147.1 

147.5 

146.7 

144.3 

142.5 

143.5 

143.3 

144.7 

141.1 

137.6 


146.7 

146.1 

145.6 

145.5 

145.2 

145.8 

145.1 

145.4 

145.8 

145.6 

145.9 

146.6 

146.9 

147.9 

148.6 

149.4 

148.8 

148.6 

148.1 

148.8 

150.1 

149.9 

151.5 

152.0 

153.5 


138.1 

139.0 

139.5 

141.7 

137.4 

134.5 

133.8 

131.9 

131.1 

131.4 

135.8 

136.4 

136.6 

138.5 

138.0 

139.4 

137.5 

136.0 

134.8 

134.3 

134.8 

130.5 

129.9 

130.4 

132.9 



10000 

OASPL 


137.3 129.0 

156.0 154.5 


135.0 136.1 152.0 

158.0 156.6 162.8 


135.6 121,5 

150.6 146.4 


Double Annular Configuration D13 (Continued) . 


713 

714 

715 

0 

Dw 

D 

Dot 

D 

D w 

140.8 

132.0 

139.8 

131.7 

138.0 

Ijx. 2 

139.3 

134.8 

139.7 

132,3 

139.3 

134.1 

139.7 

133.8 

139.2 

134.6 

139.2 

132.8 

141.9 

132,3 

140.4 

133.0 

138.7 

132.3 

139,9 

132.7 

139.1 

134.2 

137.6 

134.4 

137.0 

132.8 

138,8 

132,5 

137.0 

132.5 

134.8 

132.9 

137.1 

131.9 

136.1 

132.1 

130.6 

132.6 

135.6 

132.6 

134.9 

131.4 

130.8 

134.0 

134.6 

135.0 

134.1 

133.3 

131.4 

133.9 

136.4 

135.1 

133.4 

133.6 

133.5 

134.2 

137.5 

135,7 

135.5 

135.2 

135.4 

136.1 

138.1 

137.1 

136.6 

136.1 

135.9 

137.1 

137.6 

138.1 

137.9 

137.4 

138.0 

137.9 

139.2 

139.1 

138.2 : 

139.1 

137.2 

138.8 

140.5 

138.6 

138.0 i 

137.8 

139.4 

137.4 

140.1 

136.4 

140.4 

135.4 

137.5 

136.3 

140,2 

136.3 

139.4 

135.8 : 

135.5 

134.6 

139.2 

! 134.4 

137.7 

133.9 

138.8 

134.1 

139.3 

134.1 

138.6 

133.4 

134.8 

131.5 

138.8 

132.0 

137.5 

131.0 

134.0 

130.3 

136.5 

130.8 

134.5 

131.1 

130.5 

128.9 

134.0 

129.4 

131.0 

129.4 

127.4 

127.7 

133.4 

126.7 

128.7 

125.7 

126.1 

126.7 

■ 

125.2 

128.1 

124.2 

127.6 

126.7 

125.7 

130.4 

123.5 

132.1 

126.5 

141.3 

125.5 

135.1 

123.2 

150.9 

148.3* 

152.7 

148.6 

151.1 

148.0 







































Table VII. One-Third Octave Band Sound Pressure Levels* Double Annular Configuration D13 (Concluded) . 


RDG. 

716 

699 

723 

726 

727 

' 728 

Probe 

D 

Dw 

D 

Dw 

U 

D 

Dm 

U 

D 

Dw 

0 

D 

Dvr 

« 

D 

Dw 

frequency 

















31.5 

143.6 

131.7 

139.6 

141.0 

138.3 

136.6 

141.0 

138.3 

136.8 

138.0 

138.5 

137.6 

136.7 

138*5 

137,3 

136.5 

40 

143.3 

136.6 

139.3 

137.3 

135.4 

136.8 

138.8 

135.4 

137.0 

135. B 

134 i 9 

136.8 

135.8 

135.1 

136.5 

135.3 

50 

141,7 

135.1 

141.5 

141.6 

136.9 

138.0 

138.6 

136.4 

137.7 

137.6 

137.4 

139,5 

137.8 

137.4 

140.0 

138.6 

63 

140.2 

134.8 

145.2 

142.5 

138.1 

140.7 

138.0 

137.1 

139 .9 

138.5 

137.1 

140.9 

139.0 

137.9 

141.2 

139.0 

SO 

13B.9 

137.4 

144.6 

146.4 

140.0 

140.1 

139.4 

139.8 

139.6 

140.2 

139 . 3 

139.9 

139.4 

139.5 

140.6 

139,2 

100 

13B.3 

134.8 

142.5 

144.3 

137.4 

138,8 

139.5 

136.9 

137.8 

139.8 

137.4 

138.0 

138.5 

137.1 

138,0 

138.3 

125 

136.6 

133,6 

142.8 

142.6 

141.0 

139.8 

139.9 

139.8 

137.8 

140.9 

139.8 

138.1 

139.6 

139.5 

138.1 

138.4 

160 

134.9 

132.6 

144.1 

142.1 

138.3 

141.9 

141.1 

137,5 

140.4 

141.1 

137.5 

139.4 

140,6 

137,5 

139.1 

139.9 

200 

133.3 

135.3 

144,8 

142.3 

141.2 

140.8 

140.3 

141,5 

140.3 

141.3 

141.2 

139.8 

139.8 

141.2 

139.8 

139.3 

250 

134.2 

135.6 

146.7 

143.9 

139.3 

143.4 

142.1 

138.8 

142.7 

142,6 

138.0 

142.2 

141.1 

137.8 

140.9 

140.6 

315 

135.3 

136.7 

147.3 

145.7 

145.7 

145.3 

144.2 

144.9 

144.5 

144.2 

143.7 

142.8 

143.6 

142.9 

141.5 

141.7 

400 

134.9 

137.6 

149.6 

147.4 

152.4 

148.6 

146.1 

151.7 

148.4 

147.1 

150.2 

146.4 

144.9 

148.4 

144.4 

143.6 

500 

135.1 

138.1 

IS 2. 2 

149.4 

153.7 

149.9 

148.1 

155.2 

150,6 

150,1 

156.2 

151.1 

149.9 

151.2 

145,6 

146.6 

630 

137.7 

140,1 

154.7 

151.4 

154,0 

151.5 

150.1 

154.0 

150.2 

150.4 

154.5 

151.7 

150.6 

153.5 

149.5 

148.4 

800 

136.0 

138.8 

154.7 

151.6 

157.2 

153.5 

150.3 

154.7 

151.5 

150. E 

155.2 

150.7 

149.3 

156.7 

151.2 

150.1 

1000 

138,6 

136.4 

153.9 

152.4 

157.4 

153.1 

151.4 

157.1 

152.6 

151.9 

154.9 

150.6 

149.9 

156.1 

150.1 

148.4 

1250 

137.8 

135.8 

153.2 

149.3 

154.1 

154.0 

149.8 

155.1 

154.2 

150.3 

155.1 

153.0 

149.6 

154.4 

150,7 

147.8 

1600 

137.2 

134,4 

152.2 

149.6 

150.7 

152.7 

150.1 

151.2 

151,0 

150.6 

151.4 

152.2 

150.6 

151.4 

151,5 

149.6 

2000 

137.6 

133.6 

152.1 

149.6 

151.6 

151.1 

149.9 

152.1 

150.6 

151.4 

152,6 

150. 1 

151.4 

152.6 

149.8 

150.6 

2500 

136.0 

131.3 

151.0 

147.5 

151.3 

150.0 

146. a 

151,0 

150.5 

146.5 

150.8 

150.5 

146.8 

150.3 

151.3 

146.0 

3150 

133.5 

131.6 

149.3 

145.3 

146.9 

151.0 

146.1 

147.7 

151.8 

147.8 

147.7 

151.8 

147.8 

147,4 

152.3 

147.8 

4000 

130.8 

129.9 

148.0 

146.6 

143.3 

152.3 

146.6 

142.6 

152.5 

148.6 

142.8 

153.8 

149.4 

142.8 

153.8 

149.6 

5000 

127.4 

127.0 

148.7 

146.0 

142.8 

152.7 

145.7 

142.8 

153.4 

147.5 

142.3 

153.9 

147.0 

142.6 

154,2 

146.7 

6300 

127.1 

125.0 

144.6 

138.2 

147.5 

148.6 

137.7 

137.5 

148.9 

140.0 

137.2 

149.6 

139.7 

137.7 

149.9 

139.7 

8000 

130.4 

124.7 

142.6 

134.5 

133.3 

145.1 

133,2 

133.1 

144.9 

135.2 

132.8 

144.6 

136.5 

133.1 

143.4 

136.2 

10000 

134,6 

123.5 

142.6 

132,2 

134.3 

143.6 

132,5 

133,5 

143.8 

132.7 

133.5 

142.6 

134,2 

133.8 

142.6 

135.5 

0ASPL 

151.8 

149.3 

163.7 

161.1 

164.2 

■163.3 

160.2 

164.0 

163.0 

161.0 

163.8 

163.0 

160.5 

163.4 

162.5 

159.6 
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Figure 24. Double 



, D13, Combustor Noise Reading 706. 
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Figure 29, Double Annular, D13, Combustor Noise Reading 711. 
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Figure 44. Double Annular, D13, Combustor Noise, Upstream Probe Reading 728. 








SECTION 6.0 


DISCUSSION OF RESULTS 


Under Phase I of this program six combustors were tested and two addi- 
tional combustors have been tested under Phase II. Comparisons of the 
noise characteristics of the Phase II data with five of the six combustors 
from Phase I were performed. Effects of the various aerodynamic parameters 
were observed. There was also an additional set of data introduced by the 
downstream wall static probe utilized in the Double Annular Configuration 
D-13 test. An evaluation of this data is also presented. 

The comparison of the upstream to downstream probe spectra was illus- 
trated in Section 5.2 (Figures 15 to 22). There appeared to be good agreement 
at the lower frequencies and peak levels. However, above 1000 Hz there was 
some divergence in the levels seen by the two probes. The difference in the 
spectral shapes is believed to be due to the propagation path differences for 
each probe. The downstream probe was directly exposed to the combustor noise 
source while the upstream probe saw the blockage of the combustor dome and 
cowl and the impedance change due to an area change in the step diffuser. 

Thus, the following analysis deals primarily with the probes downstream of 
the combustor. 

In some spectra there was what appeared to be a pure tone which usually 
occurred in the 500 Hz band. Examples can be seen in Figures 18 and 25. It 
appears to be a damped resonance in the combustion system. However, it could 
nOt be ascertained whether it was caused by the combustor or the test 
facility. This analysis will not be concerned with resonances but with the 
broad band noise. 


6.1 SOUND PRESSURE LEVEL SPECTRA 


The combustor noise spectra were examined with the intent of determining 
the effect of probe type, fuel-air ratio, and configuration change on the 
noise levels. A comparison to the engine prediction spectra was also made. 

Configuration D-13 had two probes downstream of the combustor, the water 
cooled waveguide probe which was immersed in the flow and the wall probe 
which was flush with the air flow path. Spectral comparisons of the two 
probes can be seen in Figures 23 to 47. The majority of the wall data agree 
with the probe within a few dB. In Figure 46 the difference in OASPL was 3.1 
dB. However, a test condition did occur, Figure 47, where the wall acoustic 
signal appears to be overcome by the local boundary layer, resulting in a 
flatter spectrum. The levels of the two probes agree somewhat from 400 to 
2500 Hz, which provides the peak acoustic levels, but diverges at both the 
higher and lower frequencies. For readings 712, 714 and 715 at the low 3.43 
atm pressure level (Figures 32, 34 and 35) the situation seems to have 
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reversed itself. The probe levels indicate higher levels than the wall 
below 250 Hz and above 1000 Hz. With only one exception the shapes of the 
spectral data from the two probes were similar, with the wall probe overall 
sound pressure level being on the average 2.9 dB less than the water cooled 
waveguide probe. The evaluation of combustor noise was based on the signals 
from the downstream water cooled probe in order to be consistent with the 
data previously recorded. 

The effect of the fuel-air ratio can be seen in Figures 48 to 51 for 
both approach and takeoff inlet conditions at both the Phace I and Phase II 
pressure ranges. As the fuel-air ratio is increased the noise levels appear 
to increase in some regions of the spectrum and decrease in others; however, 
the net effect on the overall power level, Figure 52, indicates a consistent 
increase in noise level with an increase in fuel-air ratio for the same inlet 
conditions. At the fuel-air ratios below 0.014 the effect of this parameter 
was masked somewhat by the presence of pure tones at the 6.83 atm pressure 
and by other broad band noise sources at the 3,44 atm pressure level. It can 
never the less be concluded that the fuel-air ratio or its resulting tempera- 
ture rise should be Included into a correlation which predicts the overall 
power level. 

Figure 52 also indicates that the Double Annular D-13 takeoff levels 
are consistently higher than the CF6-50, while the approach levels appeared 
to fall on the same line. Comparing the takeoff spectra, Figure 53 (for 
the same 0.0139 fuel-air ratio) shows a greater increase in D-13 level over 
the CF6-50 occurring at frequencies greater than 630 Hz, With the fuel-air 
ratio Increased to 0.0229 (Figure 54) the difference is uniform across the 
spectrum. At this higher fuel-air ratio the spectra were roughly similar 
and peaking In the 630 to 1000 Hz range. 

Since the Double Annular combustor was tested at the lower Phase I 
(Reference 3) pressure levels also, it can be directly compared to spectra 
for the Double Annular combustor Configuration 11—11. At both takeoff and 
approach, Figures 55 and 56, there is good agreement for both the spectral 
shape and level. 

The two combustors were compared to the engine prediction spectrum and 
spectral envelope (Reference 4). The prediction spectrum had been used for 
predicting combustor noise from an engine. The actual engine data were 
found to fall in a range between the prediction spectrum and the spectral 
envelope which forms an upper bound to the range of data. The frequency at 
which the peak was set for engine data was 400 Hz ± a l/3-octave band. From 
examining the component combustor spectra it was observed that indeed the 
spectra do fall in the same shape range between the prediction spectrum and 
spectral envelope. However the frequency at which the peak occurred was 
1000 Hz, for all of the data from both the CF6-50 and D-13 (see Figure 57). 
There was good agreement from 200 Hz to 4000 Hz. The combustor noise peak 
occurring at 1000 Hz was in general true for the configurations tested under 
Phase I. Some typical examples were shown in Figure 58. 
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Figure 49. Effect of Fuel-Air Ratio on D13 at an Approach Inlet Condition. 
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Figure 58. Comparison of Phase I Configurations at Fuel-Air 
Ratio of 0.0245 to the Prediction Spectra. 
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Some of the measured combustor noise spectra could be called double 
humped. Of the two humps the major peak was the broad band hump centered 
around 1000 Hz. This was believed to be related to the combustion process 
since its level varied with fuel-air ratio. The second peak was lower in 
frequency, in the 31.5 to 200 Hz range. It was not possible to determine 
the cause of this hump. These lower frequency pressure levels were not 
believed to be combustion related since they were not influenced by the fuel- 
air ratio as the 1000 Hz hump was. The hump may be related to flow noise 
generated in the test rig. There were occasions, such as when the fuel-air 
ratios were below 0.012, that the 1000 Hz levels were as low as the lower 
frequency hump. When this occurred the overall pressure levels were in- 
fluenced to some extent. This influence was ignored in the analysis. It 
was, however, kept in mind that at the lower fuel-air ratios a noncombustion- 
related source could Influence the overall level. 

For component noise prediction purposes it is suggested that the engine 
prediction spectrum also be used for the prediction of combustor noise from 
a CF6-50 size combustor. The frequency at which the peak should be set is 
1000 Hz. The difference between the peak frequency observed in engine data 
and the component peak frequency may be attributed to the effects of the 
turbine blade rows and nozzle end impedance. 


6.2 ACOUSTIC POWER LEVELS 


In order to provide a common basis of comparison of the combustor noise 
source, it was decided to evaluate the noise on a power level basis. The 
reason for using the power level was the independence of position or location 
However, the use of the power level implied a knowledge of the source of the 
sound. It was assumed that the noise was generated in the combustor and pro- 
pagated In the plane-wave mode past the downstream probe. In an engine this 
would be the acoustic energy that would propagate through the downstream 
turbine and exhaust nozzle to be radiated to the far field as core noise. 

The acoustic power levels were calculated based on the measured downstream 
sound pressure levels which were corrected for the local Mach number, pres- 
sure, temperature and annulus area. The overall power levels are based on 
the overall sound pressure levels and corrected In the same manner. The 
details of the calculation can be found in Appendix B. 

The overall power levels of the combustors from Phase I (Reference 3) 
and the current tests are plotted as a function of the engine prediction 
parameter, (Reference 4) in Figure 59. It was observed that the component 
combustor noise levels were consistently higher than the turbofan engine data 
and they appeared to exhibit a somewhat flatter slope than the engine data. 
This suggested 1) that there was an attenuation of the combustor noise levels 
as they traveled from the combustor through the turbine blade rows, exhaust 
nozzle and jet and/or 2) the prediction parameter for component combustor 
noise is different from the parameter used to predict core noise from engines 
possibly because of the absence of the turbine, core nozzle and jet. Thus it 
may be necessary to determine first a correlation for the combustor noise at 
the source and then the effect of propagation through the attenuation mecha- 
nisms in the engine. 
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The problem of predicting combustor noise using various steady-state and 
fluctuating aerodynamic parameters has been examined in the literature (Refer- 
ences 3 and 4) . However, in order to arrive at a prediction capable of being 
used by an engine designer during the preliminary design of an engine it would 
be advantageous if only the steady-state parameters determinable from an en- 
gine cycle could be used. It was a] so initially desirable to hold some 
geometry parameters constant. As a consequence of the basic pollution pro- 
gram that required placement of the final ECCP combustor in a CF6-50 engine, 
the combustor diffuser duct and the CF6-50 combustor envelope were maintained 
constant for all of the combustors reported under Phase I (Reference 3) and 
Phase II. Therefore only the details of the combustor geometry were changed 
thus eliminating other extraneous effects. The data from Phase I involved 
some significant geometry changes. There were large variations in the way 
fuel entered the combustion 2 one and how it was mixed. Under Phase II, only 
two combustors were tested, the CF6-50 and Double Annular Configuration D-13. 
Of the two, Configuration D-13, was tested over a wider range of inlet tem- 
peratures, pressures and flow rates than any of the previous combustors. It 
xtfas felt that a correlation of the noise to the steady-state aerodynamic 
parameters based on this combustor could be used 1 ) as a basis to evaluate 
the effect of modifications in the configuration and 2 ) to provide a prelimi- 
nary design prediction parameter. 

The parameters that were selected for a final correlation with the over- 
all power level were the air flow rate, W, the temperature rise, (T 4 -T 3 ), and 
the inlet temperature, T 3 , and the inlet pressure P 3 . These variables are 
sufficiently easy to determine in an engine cycle for preliminary design 
predictions. 

A statistical correlation involving a linear multiple regression analysis 
was applied to the set of D-13 data. The result of the correlation, was a 
parameter of the form: 


OAPWL =102.7+10 log 


10 


[« 


w 1 - 29 (t 4 -t 3 ) 0 - 79 I 3 °- 48 ] 


( 1 ) 


The standard deviation from the mean line (Figure 60) was 0.90 dB. The 
correlation did not require the effect of an inlet pressure to be included 
separately. This may be because the inlet pressure and temperature are 
related through the ideal gas law even though they were varied independently. 
The air flow rate and temperature rise are each approximately to the first 
power. The product of these two terms is directly related to the heat 
release rate. 

The correlation Equation (1) from D-13 was compared to some of the other 
combustor configurations. It was observed, Figure 61, that each configuration 
fell into one of two categories. The first category (Configurations 11-11, 
11-12 and 1-16) consisted of data that fell on a single line parallel to the 
D-13 line but separated by some constant distance. The second category 
(Configurations 1-12, III-l, and CF6-50) fell on two parallel lines each 
separated by a different constant distance. Of the second category the points 
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at the approach inlet condition were closer to the D-13 correlation line. 

The takeoff inlet condition noise for those three combustors was noticeably 
lower than both the D-13 line and their corresponding approach levels. It is 
suggested that both the difference between configurations and the difference 
associated with the changed inlet conditions are related to the structure 
and level of the turbulence generated in each combustor. It is further sug- 
gested that the turbulence and its related noise from combustors of the second 
category is more sensitive to Its operating "rendition than combustors of the 
first category. 

Any general similarity of each category that was exclusive of the other 
was not immediately evident in the examination of the configurations. 

The exponents of the terras in this prediction equation were rounded off. 
Then the equation becomes: 


OAPWL = 100.2 + 10 


log M [w (T 4 -T 3 )^7] 


( 2 ) 


The standard deviation from the mean line, Figure 62, was increased to 1.2 dB. 
The trends and relationships of the other configurations remained essentially 
unchanged. The simplicity of this equation was attractive and it is recom- 
mended as a preliminary design prediction. It should be noted that this 
equation was based on component combustor data. The effect of a change in 
end impedance from that which existed in the component test to that which 
occurs in an engine is not included. 
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SECTION 7.0 


ATTENUATION THROUGH TURBINE BLADE ROWS 


The attenuation of combustor-generated noise by an engine turbine was 
suggested in Reference 4. An analysis of the attenuation phenomenon was 
performed by Bekofske (Reference 5) for subsonic relative flow. This analysis 
is extended here to include the case where there is subsonic relative flow 
at the inlet of a blade row but supersonic at the discharge! side. The 
turbine attenuation for the CF6-50 engine was calculated based on this 
analysis. The calculated attenuation was compared to an estimated value of 
turbine attenuation. The estimated value was arrived at by first estimating 
the power level of combustor noise in the engine based on Equation 2. The 
CF6-50 engine core noise power level based on far- field measurements was- 
then subtracted from the component level. The resulting difference was 
called the turbine attenuation and was compared to the analytic solution. 

The comparison showed that the analytic solution over predicts the 
attenuation. 

In addition the analytic prediction method was exercised in a parametric 
study and a turbine pressure ratio dependency discerned for the plane-wave 
case. 


7.1 EXTENDED ANALYSIS 


Introduction 


The attenuation of acoustic energy of low frequency waves propagating 
through a gas turbine stage was analyzed. An earlier investigation of this 
blade transmission problem (Reference 5) assumed that the relative Mach 
number was subsonic at both the inlet and discharge sides of the nozzle and 
bucket row. Results indicated that the attenuation of acoustic energy 
across typical turbine stages varies from less than 1 dB to more than 10. 

The present study extends the previously developed theory to account for 
supersonic relative flow at the discharge of a blade row. Such a choked 
condition occurs typically In single-stage, high-pressure turbines and 
highly loaded fan turbines at high power settings. The theory is derived 
for an isolated blade row. As in the earlier study (Reference 5), Mach 
number changes and flow turning across a stage element are taken into 
account. 

In actual practice, most of the energy-containing disturbances arising 
in the combustion chamber of a gas turbine have a characteristic acoustic 
wavelength which is large compared to both the chord and transverse spacing 
of a blade row. Core engine noise spectra are presented by Emmerling 
(Reference 4) and Pickett (Reference 6), which indicate a low frequency 
peak, identified as combustor noise, centered at about the 400 Hz 1/3-octave 
band. It is suitable to investigate the transmission of such low frequency 
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combustor disturbances through a turbine stage element (nozzle or bucket 
row) using a two-dimensional (cascade plane) actuator disk model. Therefore, 
as in the previously mentioned /study (Reference 5) , the analysis will be 
restricted to low frequency disturbances and the actuator disk model will be 
employed. 

One of the earlier uses of the actuator disk model was by Rannie and 
Marble (Reference 7) . Since the fluid motion in axial flow turbomachinery 
is three-dimensional, and so complex, it is necessary to use some simple 
mathematical model in most practical irstances. The actuator disk model to 
be employed here uses a concept borrowed from propeller theory and applied 
to a blade row within an annular passage. It is assumed that the radial 
flow is negligible, thus reducing the problem to one of. two dimensions (the 
cascade plane). Further, the actuator disk is a model in which each blade 
row is assumed to be an appropriate plane of discontinuity in the flow. 

Thus, the actuator disk method is restricted to low frequency disturbances 
whose wavelengths are large compared to the chord and transverse spacing of 
the cascade blade row. A comprehensive discussion of the actuator disk 
theory in turbomachinery is given by Hawthorne and Ringrose (Reference 8) , 
Pickett (Reference 6) employs the actuator disk model to study the inter- 
action of vorticity and entropy disturbances arising in the combustion 
chamber with a blade row. Cumpsty and Marble (Reference 9) also use the 
actuator disk model to investigate the generation of noise by entropy fluctua- 
tions into a turbine. 

The problem formulation follows closely an unpublished theory by 
Dr. R. Mani for the discharge reflection coefficient from a blade row at low 
frequencies and the analysis contained in Reference 5. A sound wave which 
impinges on the upstream side of a blade row will lead to a reflected wave 
on the upstream side and a transmitted wave and a vorticity wave on the 
downstream side of the blade row. The equations necessary to solve for the 
unknown wave-amplitude coefficients are obtained by requiring continuity of 
mass flow and energy and, depending on whether the Mach number is subsonic 
or supersonic at the discharge of the blade row, either applying the Kutta 
condition at the trailing edge of the blade row or applying the choking 
condition to the blade row. 


Problem Formulation 

The actuator disk model is shown in a two-dimensional flow field in 
Figure 63. A Cartesian coordinate system is fixed with respect to the blade 
row (nozzle or bucket) . The x-axis is in the axial flow direction and the 
y-axis is tangent to the blade row and positive in the direction of bucket- 
rotation. It is assumed that the axial Mach number is subsonic on both 
sides of the blade row. Also, it is assumed that the inlet relative Mach 
number is subsonic whereas the exit relative Mach number can be either 
subsonic (Case I) or supersonic (Case II). All quantities (angles and 
velocities) are taken to be positive when in the direction indicated in 
Figure 63. The angles a and g are the incident and exit relative stream 
angles, respectively, measured from the x-axis. Upstream and downstream 
flow parameters will be indicated by subscripts 1 and 2, respectively. 


00 

00 

Relative Mach Numbers Indicated: 

M l» M 1X> M 2X ^ 1 . 

Case I Mg <C l 

Case II Jl 2 )1 



Reflected 



Figure 63. Geometry of Wave 


(Direction of Bucket Rotation) 
v 



Incident on Stage Element. 



Wave Description 


Consider a sound wave with frequency, to, and wave number, k]_, which 
impinges on the upstream side of the blade row at an incident angle, 0j. 

Such an incident wave will lead to a reflected wave at angle 0p v and a trans- 
mitted wave at. angle 9y. The angles 0p and 0^ are to be calculated. In 
addition, a vorticity wave is induced on the downstream side of the blade 
row. 

The physical mechanism by which a "wake" of vorticity is produced by an 
airfoil in no nuniform motion is described by Karman and Sears (10) . Briefly, 
the unsteadiness in lift gives rise to a wake composed of continuously dis- 
tributed vortex lines behind the airfoil. Within the approximation of the 
actuator disk theory (infinitesimal blade spacing) a coalescence of vortex 
wakes (vorticity wave) is shed from the trailing edge of a blade row on 
which unsteady distrubances have impinged. 


The frequency of each of the above mentioned waves is to. However, 
because the temperatures (and, therefore, the speeds of sound) are different 
upstream and downstream of the blade row, the wave numbers will he different. 
The wave number is related to the frequency by k = co/a where a is the ambient 
speed of sound. Therefore, the upstream and downstream wave numbers are 
related by 


k 2 aj 
ki W2 


(3) 


The pressure disturbance associated with the incident acoustic wave can 
he written as 

k^(x cos Qj -H y sin 0^.) 

Pi = p ex P j h + H, cos 07 + "5; ' sto'e" * “ t} C4: 

.... lx I ly 1 

where P is the amplitude of the incident wave, M]_ x ana M]_y are axial and 
tangential components of the upstream relative Mach number, and j = J-L. The 
prime indicates a small deviation from the undisturbed flow field. 

The pressures associated with the reflected and transmitted waves can be 
expressed as 


p^ » RP exp j {yr' 


k 1 (— x cos 0 ^ + y sin 0 R ) 


R 


and 


^lx COS ^R + M iy 8111 


- cot) 


k 2 (x cos 0 T + y sin 0^,) 

p' = TP exp j — rrr — 5 — — - - tot} 

T r J 1 + cos Q_ + M_ sin 

2x T 2y T 


where R and T are unknown coefficients. 


(5) 


( 6 ) 
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The tangential or y dependence of the vorticity wave must be the same 
as that of the incident wave so that it has the form 


exp j {k^x + k^y - mt} 


k^ sin 0-j. 


1 + M ix cos 0 x + M iy sin 9l 


The x dependence of the vorticity wave is determined by the fact that, in a 
frame of reference fixed with the fluid (xp = x ~ a£ M 2 x t; yp = y - a2 M 2 yt) , 
such a wave appears stationary; that is, its time dependence vanishes* In 
this frame of reference moving with the fluid, expression (7) is written 


exp j{k x (xp + M 2 x^ + k y^ y P + a 2 M 2y^ " 


Since the time dependence must vanish, 


k x a 2 M 2 x + k y a 2 U 2y ~ “ “ 0 


k 2 k y M 2 v 


The vorticity wave is divergence free, so that the velocity components of 
the vorticity wave must satisfy 

3u* 3v’ 

+ = 0 . (113 

3x 9y 


Using the expression (7) in equation (11) , it follows that 

v* k 

= _ 

u* k 

Q y 

Therefore, the velocity components of the vorticity wave are given by 

(u i- V = exp 3tk * x + V ' “ t} 

x y x y 




where Q is an unknown amplitude coefficient and ky and k x are given by 
equations (8) and (10). There are no pressure or density fluctuations 
associated with the vorticity wave, so that 

p' = o 

Q U 


P Q = ° 


Velocity and density fluctuations associated with the incident, reflected 
and transmitted waves are given by 

(u r V P = ( cos e r sl " V 5^ 

k^(x cos + y sin 0^.) 

exp + m cos e -p m sin 0 T 
lx I ly I 


Pj “ exp j{- 


k^(x cos 0^. + y sin 8^) 
r+ CDS 0J. + 1L sin 0J. 


- tot}. 


(u^» V^) - (“COS 0 R , sin 0 R ) 


P 1 a l 


exp j{- 


k- (-X cos G n + y sin 0„) 
J. ft ft 


, cos + M, sin 0^ 
Lx R ly R 


, RP , , 1 

p^.-jaxp j{— 


k n (-X cos 6- + y sin 0._) 
J. ft ft 

“ M- cos 0,, + sin 0, 
lx R ly J 


- tot } , 


“ tot}. 


(ujL, v') “ (cos 0 , sin e ) - — — 
T I i T P 2 a 2 


exp j{ 


k^(x cos 0 T + y sin 0 T > 

1 + cos 0 ' + M 0 sin 0„, 
2x T 2y T 


“ tot} 


Tp k 2 (x cos 0 T + y sin 0 T ) 

p i = 4 exp J{ i + a 2x co S e T + M 2y si„ 8 t 


“ U)t}. 




Relations expressing 0 R and 0^ in terms of 0y are found using the fact 
that the y dependence of the reflected and transmitted waves must also coincide 
with that of the incident wave. Therefore, from equations (4) - (6) 
sin 0^ k 2 sin 0^ 

1 - «XX cos 0 E + 'Ay Sin e R “ 1 + M 2x cos S T + M,, y sin 0 T 


lc^ sin 0-j. 

1 + M, cos 0 T + M, sin 0 T ^2 
lx I ly I 


(19) 


After some manipulation, it follows that 

1 (1 - M? ) sin 0 T 

0 R = tan [- ~ — ' ] 

<1 + M^)cos e i + 2^ 


and 


0 T = tan 


-1 ®2x (1 - '"V + G ^ (1 " ®2y>" - Gi(1 - 


2y' 

(1 - GM 2y > 2 ” G 2 


( 20 ) 


( 21 ) 


where G(0j) is given by equation (19). The reflected and transmitted waves 
should always appear on the opposite side of the machine axis from the incident 
wave. Therefore, with the notation and sign conventions of Figure 63, if 
0-1 > 0, then Qr > 0 and G-j > 0, and if 0j < 0, then 0 r < 0 and 6j> < 0. If 
0£ = 0, then 6 r ~ = 0. 


Only the acoustic waves for which the normal to the wave crest satisfies 
j 0 I l < 90° + sin" 1 M lx (22) 


impinge on the blade row. Other waves are propagating upstream faster than 
they are convected downstream by the flow. Similary. only those transmitted 
waves for which 

| 0 T | < 90° + sin” 1 M 2x (23) 

are propagating downstream. 

When the radical in equation (21) vanishes, i.e., 

(1 - GM 2y ) 2 - G 2 (l - M 2x ) « 0, (24) 
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( 25 ) 


equation ( 21 ) leads to 
tan 9 T, cut-off = 


-M 


' 2 x 


where the sign in the denominator has been preserved to identify the quadrant 
of Of. Equation (25) expresses the same limit on ©t as given in the inequality 
(23). When the radical in equation (21) becomes negative, the transmitted 
wave given by equation ( 6 ) is decaying exponentially. This downstream cut- 
off phenomenon leads to a further restriction on the incident angle, 0 i. 

Solving equation (17) for 0j = fen ( ©'j* ) and substituting equation (25) gives 
01 cut-off* Acoustic waves whose angles are beyond ©X, cut-off can impinge 
on’ the blade row but are essentially shielded because of condition (23) . 


Governing Equations 

Three equations involving the disturbed flow parameters are necessary 
to find the unknown amplitude coefficients for the generated waves, T, R, 
and Q. The flow disturbances must always conserve mass and energy which 
provide two matching conditions, furthermore, for the case of subsonic 
relative velocities on both sides of th blade row, the Kutta condition can 
be applied at the trailing edge of the blade row (Case I). For the case of 
subsonic relative velocity at the inlet but supersonic relative velocity at 
the discharge side of the blade row, the choking condition can be applied to 
provide the third equation (Case II) . 

The linearized equation for the continuity of mass flow is given by 


[Up 1 .+ pu' ] 1 = [Up 1 + pu T ] 2 


(26) 


The subscripted square brackets indicate evaluation of the quantity on the 
upstream (x — 0 “) and downstream (x = 0 "*") sides of the blade row, respectively. 
The linearized equation for the conservation of energy per unit mass, when 
written in a frame of reference fixed to the blade row, can be written 

fp— + Uu’ 4- Vv 1 ]^ = + Uu* + Vv’ j 2 -. (27) 

Case I, subsonic relative velocities throughout: the Kutta condition 

requires that at the trailing edge the flow should be tangent to the exit 
relative stream angle, 0 (see Figure 63) . The unit vector normal to the 
flow is 


e - - sin (3 e + cos 3 e . 
n x y 

Then the Kutta condition requires that, at the trailing edge. 



(28) 


(29) 
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Case II, supersonic relative velocity at the blade row exit: this case 

assumes that no shock waves occur in the blade passage; the presence of shock 
waves would allow a subsonic relative velocity at the discharge. The choked 
condition will manifest itself in the fact that the corrected mass flow is 
dependent only on the upstream stagnation parameters dD , A discussion of its 
application to disturbed flow at a blade row can be found in Reference 9. The 
choked condition requires that 

= Constant 



where W = mass rate of flow, 

A — cross sectional area, 

R = gas constant, 

Y = ratio of specific heats. 


Tq =' stagnation temperature, 
and Pq = stagnation pressure. 


The mass rate of flow is w = pAU, For y = constant, which will be assumed in 
this analysis, Equation (30) can be Witten 


oIJ 



constant. 


(31) 


Condition (31) is used to obtain the third equation involving disturbed 
flow parameters when the relative velocity at the blade row exit is supersonic. 
Introducing 




into Equation (31), where is the absolute Mach number, and taking the 
logarithmic differential yields 



(32) 




where 


T* - temperature disturbance, 


and 


pn + 


3L=-i M 2 


ab 


The absolute Mach number is 


M , = 
ab 




where C is the absolute flow velocity. Denoting the absolute flow velocity 
disturbance by c’, it follows that 

_ v ± 1 ML = (Y + 1)M ab ^TL _ 2 c^ 


y - 1 M 


2[l + ^ M 2 b ] 


When this expression is substituted into Equation (32), it becomes 


(Mi - 1) £- + 2K - (y + DM 2 . = 


(33) 


The temperature disturbance in isentropic flow is related to the pressure 
disturbance by 

lL = V ~ 1 p 1 

T TP* 

The absolute flow velocity disturbance, c*, can be expressed in terms of the 
relative velocity component disturbances as 

c* U* V* 

— coscj) + sin<f> 


where tf) is the absolute flow angle, measured from the axial direction. There- 
fore, equation (33) is written 

(y +• Dii 


. (r^D . M | b _ x) p . + [•§■“- - 

Cy + DM ab _ 


ab 


cos^lu 1 


(34) 


- [— 


sin<j>]v' - 0. 


It should be noted that in Equation (34), u 1 and v f can contain, in general, 
contributions due to both pressure and vorticity waves. This equation will 
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be applied upstream of the blade row and used as the third condition (choking 
condition) for the case of supersonic relative flow at the discharge of the 
blade rox*. 


When equations (4) to ( 6 ) and (13) to (18) are substituted into equations 
(26) and (27) the following two equations result: 

— [M,„ + CDS 0 t ]T + ^ 7^=5 + Icos 0 R - - + COS 0j- (35) 


a 2 ' 2 x 


l 2 

v x x 


and 


'i P 1 ' ' Q 

~ [1 + M 2 COS (3 - 0 t )]T + ~ Iky ^2x - M 2y 3 / 2 "~ — j 
2 ^X + y- 

- [1 - M cos (a + e R ) ]R = 1 + \ cos (a - Bj) 


(36) 


A third equation is obtained by substituting Equations (4) to ( 6 ) and (13) 
to (18) into either Equation (29) or (34) , depending on whether the discharge 
relative velocity is subsonic or supersonic: 

[sin(0 T - 0)]T - sinS + cos 8 ] 


x ' /FT? 

x y 


0 (Case I) 


(37a) 


or 


2 p n 


[(v + 1)M cos(^+ 0 R ) - cos 0 R + (Y - DC M aab “ *■> ^ 


lx 

2y 


= (y + 1 )M lah cosC^ - B z ) - ~ cos 0 ] . - (y - D “ D (Case II) (37b) 

lx 


Equations (35), (36), and either (37a) or (37b) can be solved for T, Q, and R. 


Attenuation of Acoustic Energy 


After the transmission coefficient, T, has been determined, it is 
necessary to relate it to the attenuation of acoustic energy. Using 
Blokhintsev's ( 12 ) results, the acoustic energy density can be written 

,2 

e = ~ P — - . 7 (a + C * e ) 
pa J ■ p . 


(38) 


where C is the absolute flow velocity and e is the unit vector normal to 
the wave front. In Equation (38), note that the pressure fluctuation down- 
stream of a blade row is related to the pressure fluctuation of an incident 
wave by the transmission coefficient, T, in Equation ( 6 ). 
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( 39 ) 


The intensity flux vector has the expression 

,2 

1 ~ "y (a + C * e ) (a e + C) « 

3 - p p 

pa 


Then the attenuation of acoustic energy across a blade row is given by 

I, . 6 

Attenuation = 10 log^Q Cj — — ') 

-1 x 


(40) 


where only the axial component of the intensity flux vector is of interest. 
The attenuation given by (40) is a function of the incidence angle of the 
wave, 0i, as well as the other flow parameters. 

A listing of the computer program incorporating this extended analysis 
can be found in Appendix D. 

A blade row attenuation computer study was conducted for a single stage 
high pressure turbine with two different inlet temperatures = 950° F and 
350° F. Subsoni c cas es (Pressure Ratio Fr = 1.9, 2.0, 2.2 and 2.5) for 70, 
90, 100, 110% N/ were run. The attenuations for normal incidence (0j = 
0) waves were used. Figure 64 illustrates the results. The overall stage 
attenuation is indicated as 8 to 9 dB over the entire pressure ratio range. 

The prediction program was then exercised for six different turbines 
over a larger range of operating conditions. There were some definite 
trends for the plane wave (zero incidence) case. The attenuation increased 
with the stage pressure ratio (total-to-static) as shown in Figure 65. The 
trend suggests using (AT) or (AT/T) as the correlating parameter where: 

Y“1 



and T is the turbine inlet temperature. 

The last stage of a low pressure turbine offers significantly less 

attenuation than the preceding stages for the same pressure ratio. This 

would suggest that the blade turning or exit angle plays a role, but the 

role is not immediately obvious. The analysis suggests a similar result 

when a plane wave assumption is made. Using 0j = ... 0, 9r = ®T = ky = 0 and 

(k x / /k.2-rk^) = 1, it can be shown that 
x y 













and 




(42) 


where 

T = (transmitted/incident) acoustic pressure 

A = attenuation in dB 

—1 

8 = exit swirl (tan“ Mo /M„-) 

y 

The quantities (a-^/^) and CP 1 /P 2 ) are related to the work extraction by the 
turbine and hence to the pressure ratio, Pr. Furthermore, since the axial 
velocity through a blade row is maintained approximately constant, the ratio 
t(l + M^ x )/(1 4- M 2 x )] also becomes, roughly speaking, a function of Pr! 
Equation (42) can then be recast as 

A = function (Pg, 3) 

The current prediction method for aircraft engine core noise uses 
three separate prediction lines for turbojets, turboshafts , and turbofans. 
The difference between the lines has been attributed, at least in part, 
to differing low frequency noise attenuation by the turbine stages. The 
pressure ratio (Pr) dependency indicated by the analysis would suggest 
inclusion of Pr or (AT/T) in the core noise prediction parameter - in 
order to collapse the three lines. 


7.2 EFFECT OF BLADE ROW ATTENUATION IN THE COMBUSTOR NOISE PREDICTION 

The CF6-50 blade row attenuation was calculated using the extended 
analysis for three engine conditions. The individual stage and total 
attenuation for the approach, cruise and takeoff conditions are listed in 
Table VIII. For this operating range, the total .calculated attenuation 
varied between 34 and 37 dB. 

For comparison purposes, the component combustor noise was calculated 
for the same three engine cycle conditions based on equation (2) . These 
extrapolated component combustor overall power levels are listed in Table 
IX. The CF6-50 engine core noise was then calculated from farfield test 
data. The engine power levels were based on fitting the combustor noise 
prediction spectrum with the peak frequency at 500 Hz to the peak angle 
(120° from the inlet) spectra. The directivity that was used was the one 
recommended in Reference 4. Subtracting the farfield level from the 
predicted component combustor level was the attenuation. At approach the 
attenuation was 14.7 dB while at cruise and takeoff it was greater than 
16.6 and 14.8 dB respectively. At the cruise and takeoff condition jet 
noise dominated the spectra. The maximum level that the combustor noise 
could assume without effecting the farfield spectra (10 dB below the 
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Table VIII. Calculated CF6-50 Turbine Stage 

Attenuation in the Axial Direction. 


Turbine 

Stage 

Approach 

Cruise 

Takeoff 

High Pressure 

1 

10.6 

11.1 

11.1 


2 

6.3 

7.3 

7.3 

Low Pressure 

1 

5.9 

5.6 

5.6 


2 

5.3 

5.6 

5.6 


3 

4.1 

4.9 

4.9 


4 

1.8 

2.4 

2.3 

Total Attenuation 

■ s; - 

34.0 

36.9 

36.8 


Table IX. Combustor Noise Attenuation in the CF6-50 Engine. 



Approach 

Cruise 

Takeoff 

Component Combustor OAPWL 

159.4 

162.8 

164.0 

Farfield Engine Core OAPWL 

144.7 

<146.2 

<149.2 

Attenuation 

14.7 

>16.6 

>14.8 



measured level) were listed. The core noise must be at that level or 
less for those two cases . 

The discrepancy between analysis and indicated experimental results 
could be due to either the unique end impedance condition involved in the 
component tests or the simplifying assumptions involved in the analysis, 
perhaps both. The latter invokes an actuator disk assumption and also 
uses an isolated stage - thus neglecting the effect of reflections from 
upstream and downstream stages, A study is currently being performed 
under NASA Contract (No. NAS3-19435) in an attempt to resolve some of 
these questions. 



SECTION 8.0 


CONCLUSIONS 


Although it remains to he proven that the pressure levels recorded dur- 
ing this test were entirely acoustic pressures, conclusions may be drawn 
assuming that the peak levels are indeed combustor noise. 

# A wall static probe may be used to measure combustor noise, 
however locally turbulent test conditions may occur which 
may generate pressures in excess of the acoustic pressure 
levels . 

£ The frequency of the peak noise was observed to be 1000 Hz 

in the component tests rather than 400 Hz ± one-third octave 
band observed from engine data (see Reference 4). 

9 The overall power level was found to be a function of the 
heat release rate and the inlet temperature. 

# The attenuation of combustor noise was found to be signi- 
ficant (greater than 14 dB). 

$ The blade row attenuation analysis as it is presently 

structured over-predicts the apparent attenuation based on 
component combustor tests . This may be due to an initial 
assumption in the analysis and/or the change in combustor 
exit impedance due to the addition of a turbine when 
going from component data to engine data. 
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APPENDIX A 


NOMENCLATURE 

a 

A 

B 

c* 

C 

e 

D 

f/a 

G 

H 

I 

j 

k 

L 

M 

M ab 

P R 

P 

pV 

P, Q, R, T 
R 

RDG 


SYMBOLS 

Speed of sound 

Cross sectional area 

Fuel nozzle spacing 

Absolute flow velocity disturbance 

Absolute flow velocity 

Unit vector 

Downstream, Engine Station 3.9 

Fuel Air Ratio 

Defined by Equation (17) 

Height 

Intensity flux vector 
vCi 

Wave number 
Length 

Relative Mach number 
Absolute Mach number 
Turbine stage pressure ratio 
Pressure (atmospheres) 

Pressure disturbance 
Amplitude coefficients 
Gas constant 
Reading number of test 


NOMENCLATURE (Continued) 

T Temperature (° 1C) 

T 1 Temperature disturbance 

t Time 

u 1 Axial velocity component disturbance 

U Upstream, Engine Station 3.0 or mean axial velocity component 

v T Tangential velocity component disturbance 

V Mean tangential velocity component 

W Mass rate of flow 

x Axial Cartesian coordinate (fixed to blade row) 

y Tangential Cartesian coordinate (fixed to blade row) 

a Incident relative stream angle 

3 Exit relative stream angle 

Y Ratio of specific heats 

A Increase or decrease 

e Acoustic energy density 

0 Wave propagation angle 

* 1+2 r 4, 

p Mean density 

p' Density disturbance 

tj) Absolute flow angle (relative to axial direction) 

m Frequency 

SUBSCRIPTS 

B Burner 

C Combustor 

D Dome 
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SUBSCRIPTS (Continued) 


F Fixed with, respect to fluid 

f Fuel 

I Incident wave 

n Normal to flow 

p Normal to wave front 

Q Vorticity wave 

R Reflected wave 

S Static 

T Transmitted wave or total 

W Wall 

x Axial 

y Tangential 

0 Stagnation value or ambient condition 

1 Upstream of blade row 

2 Downstream of blade row 

3.0 Engine station upstream of the combustor 

3.9 Engine station downstream of the combustor 



APPENDIX B 


POWER LEVEL CALCULATION 


The acoustic power levels were calculated using the measured sound 
pressure levels from a probe located in an annular duct with mean flow. 

The acoustic power levels were calculated (Reference 13) for the downstream 
probe only. For this probe the acoustic waves were traveling in the direc- 
tion of the flow and all of the energy was assumed to be concentrated in the 
plane-wave mode. If A was the annulus area in square meters at the probe 
sensing plane, the power level correction equation in dB re, 10“13 watts was: 


FWL = SPL + 9,9 + 101og l0 



+ lOlog 


10 


T 

T 0| 


+ 101og l0 [A] 


The pressures and temperatures were static but, because the Mach number was 
loxf, the total pressure and total temperature may be used without appreciable 
error. 

The above equation was derived from the definition of power. 

2 

PA 

Watts = 

PC 


The pressure squared term was contained in the SPL and the characteristic 
impedance pC was corrected to standard day conditions by: 



The correction for the flow incorporated the square of one plus the Mach 
number. This was a simplification of the correction suggested by 
Blokhintsev (Reference 12) : 

(1 + cos 0) (cos 0 4- M ^) 

For the plane-wave mode these corrections are identical. 



APPENDIX C 


NARROWBAND SPECTRA 


All of the spectra in this section are from 0 to 5000 Hz with constant 
10 Hz band width. The sample length was 12,8 seconds. These spectra are not 
corrected for probe losses. Those corrections can be found in Section 3,3. 


109 



SPL, dB re: 0.0002 



Figure C-l, CF6-50 Combustor Reading 317 
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Figure C- 3. CF6-50 Combustor Reading 318 Upstream Probe Narrowband Spectrum. 
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31S Downstream Probe Narrowband Spectrum. 



2 


3 


4 


5 


Frequency, kHz 

Figure C-5. CF6-50 Combustor Reading 3X9 Upstream Probe Narrowband Spectrum. 
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Probe Narrowband Spectrum 
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Figure C-7. CF6-50 Combustor Heading 320 Upstream Probe Narrowband Spectrum. 
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Figure C-8. CF6-50 Combustor Reading 320 Downstream Probe Narrowband Spectrum 
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Figure C- 9. CF6-50 Combustor Reading 321 Upstream Probe Narrowband Spectrum 
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Figure C-lli CF6-50 Combustor Reading 322 Upstream Probe Narrowband Spectrum. 
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Figure C-12, CF6-50 Combustor Reading 322 Downstream Probe Narrowband Spectrum 
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Frequency, kHz 

Figure C-13. CF6-50 Combustor Reading 323 Upstream Probe Narrowband Spectrum. 
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Figure C-14. CF6-50 Combustor Reading 323 Downstream Probe Narrowband Spectrum 
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Figure C-15, CF6-50 Combustoi’ Reading 324 Upstream Probe Narrowband Spectrum, 
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Figure C-I6. CF6-50 Combustor* Reading 324 Downstream Probe Narrowband Spectrum 
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Figure C-17. Configuration D13 Reading 705 Downstream Probe Narrowband Spectrum. 


Frequency, kHz 


Figure C-1S. Configuration D13 Reading 705 Downstream Wall Narrowband Spectrum. 
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Figure C- 20. Configuration D13 Reading 706 Downstream Wall Narrowband Spectrum 
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Figure 021. Conf iguration D13 Reading 707 Downstream Probe Narrowband Spectrum. 
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Figure C-22, Configuration D13 Reading 707 Downstream Wall Narrowband Spectrum. 
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Figure 023. Configuration D13 Reading 70S Downstream Probe Narrowband. Spectrum, 
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Figure C-24. Configuration D13 Reading 708 Downstream Wall Narrowband Spectrum 
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Figure C-25, Configuration D13 Reading 709 Downstream Probe Narrowband Spectrum. 
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Figure C-26. Configuration D13 Reading 709 Downstream Wall Narrowband Spectrum 
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Frequency, kHz 

Figure C-27. Configuration D13 Reading 710 Downstream Probe Narrowband Spectrum. 
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Figure C-28. Configuration D13 Reading 710 Downstream Wall Narrowband Spectrum 
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Figure C-30. Configuration D13 Reading 711 Downstream Wall Narrowband Spectrum 
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Figure C-33. Configuration D13 Reading 748 Downstream Wall Narrowband Spectrum. 
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Figure C-34, Configuration D13 Reading 712 Downstream Probe Narrowband Spectrum 
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Figure C-35. Configuration D13 Reading 712 Downstream Wall Narrowband Spectrum, 




SPL, dB re: 0.0002 D 





SPL, d B re: 0.0002 D 



2 3 

Frequency, kHz 


Figure C-37. Configuration D13 Reading 713 Downstream Wall Narrowband Spectrum. 
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Figure C-39. Configuration D13 Reading 714 Downstream Wall Narrowband Spectrum. 
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Figure C-43. Configuration D13 Readi 
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Figure C-45. Configuration D13 Reading G99 D-jwn stream Wall Narrowband Spectrum. 
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Figure C-17. Configuration D13 Heading 723 Downstream Wall Narrowband Spectrum. 
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Figure C-4S. Configuration D13 Reading 726 Downstream Probe Narrowband Spectrum. 
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Figure C-49, Configuration D13 Reading 726 Downstream Wall Narrowband Spectrum, 
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Figure C-50, Configuration D13 Reading 727 Downstream Probe Narrowband Spectrum. 
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Figure C-51. Configuration D13 Reading 727 Downstream Wall Narrowband Spectrum. 
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Figure C-52. Configuration D13 Reading 728 Downstream Probe Narrowband Spectrum. 
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Figure C~53. Configuration 1)13 Reading 728 Downstream Wall. Narrowband Spectrum. 



APPENDIX D 


BLADE ROW ATTENUATION COMPUTER PROGRAM 


The computer program listed here is based on the extended analysis 
described in Section 7.1. 



OOOIO REAL K1 , K2,KRA,KX0VKY, Ml t M2,MI A, M2A,,MIT,M2T 
02Q REAL M3 A, M4A, M3T, M4T 

j ) DIMENSION QQ( 100) ,A1TEN1 ( 100) , ATTEN2 Cl 00), QGSAVE( I OO) 

00031 DIMENSION XIC0RD( 1 00) , X2C0RD< I 00) ,XI (lOO) t X2( 100 ) 

00032 DIMENSION X3( 1 00) ,Y1 ( 1 00) ,Y2( 1 00) , Y3U GO) 

00035 INTEGER E,F 

00040 PI=3. 1415927 5 GR AV=32 . 1 7405 ? C0NV=778 . 26 
00050 GAMMA* 1 .32 !RGC=0. 06854 

00060 PR INI * 11 ABSOLUTE VELOCITIES ARE INPUT IN OPTION 1 " 

00070 PRINT* 11 RELATIVE MACH NUMBERS ARE INPUT IN OPTION 2” 

00080 PR I NT* “SELECT OPTION 1 OR 2 " SREAD* INPUT 

00090 IF( INPUT .HQ. 2) GO TO 130 

00100 PRINT*" ABSOLUTE VELOCITIES ARE IN FEET/SECOND" 5 PR I NT 1460 
00110 PRINT*" ABSOLUTE FLOW ANGLES ARE IN DEGREES RELATIVE TO FACE PLANE" 
00120 PRINT*" AND NEGATIVE IF FLOW DIRECTION OPPOSITE WHEEL ROTATION" 
00130 130 PRINT 1460 » PRINT* "PRESSURES ARE IN PSIA" 

00140 PRINT 1460 SPR INT* "TEMPERATURES ARE IN DEGREES RANKING" 

00150 150 N=0 SPRINT 1460 
00} 60 DO 170 1=1,100, I 
00170 170 ATTEN2 ( I ) =0 . 

00180 IFC INPUT .EQ. 2) GO TO 290 

00190 PRINT*" INPUT DATA UPSTREAM OF NOZZLE" 

00200 PRINT* "AXIAL VELOCITY, ABSOLUTE FLOW ANGLE" SREAD* V1A.PSI1 
00210 PRINT* "STATIC PRESSURE, TOTAL. TEMPERATURE" SREAD* PS! ,TZ1 
00220 PR I NT.* "INPUT DATA BETWEEN NOZZLE AND BUCKET" 

00230 PR INT* "AXIAL VELOCITY, ABSOLUTE FLOW ANGLE" SREAD* V2A,PSI2 
'0240 PR INT* "STATIC PRESSURE, TOTAL TEMPERATURE" {READ* PS2,TZ2 
.0250 PRINT* "INPUT DATA DOWNSTREAM OF BUCKET" 

00260 PRINT*" AXIAL VELOCITY, ABSOLUTE FLOW ANGLE" SREAD* V3A, PSI3 
00270 PR INT* "STATIC PRESSURE, TOTAL TEMPERATURE" ;READ*PS3, TZ3 
00260 PRINT*" INPUT WHEEL SPEED(FPS)" SREAD* UW23 SUW=0. 

00290 29Q IF (INPUT .EO. I) GO TO 420 
00300 PR INT*" INPUT DATA UPSTREAM OF NOZZLE" 

00310 PRINT*" AXIAL AND TANGENTIAL RELATIVE MACH NUMBERS'* SREAD *M1 A, MIT 
00320 PR I NT * "STATIC TEMPER ATURE AND PRESSURE" SREAD* TEMPI , PS 1 
00330 PR INT*" INPUT DATA DOWNSTREAM OF NOZZLE" 

00340 PRINT *" AXIAL AND TANGENTIAL . RELATIVE -MACH, NUMBERS" ?READ*M2A, M2T 
00350 PRINT* "STATIC TEMPERATURE AND PRESSURE" SREAD* TEMP2,PS2 
00360 PRINT*" INPUT DAT i. UPSTREAM OF BUCKET" 

00370 PRINT* "AXIAL AND TANGENTIAL RELATIVE MACH NUMBERS" SREAD* M3 A, M3T 
00380 TEMP3*TEMP2 SPS3=PS2 

00390 PRINT* "INPUT DATA' DOWNSTREAM OF BUCKET" , " 

00400 PR I IDT*" AXIAL AND TANGENTIAL RELATIVE MACH NUMBERS" S READ * M4A , M4T 
00410 PR I NT* "STATIC TEMPERATURE AND PRESSURE" SREAD* TEMP4,PS4 
00420 420 IF (INPUT .EQ. 2) GO TO 440 . 

00430 PS 1 1 =PI*PSI 1/180. SPSI2*PI*PSI2/180. ?PSI3=PI*PSI.3/l 30. 

00440 440 N-N+1 SlF( INPUT .EQ. 2) GO TO 530 
00450 V 1 T=V 1 A/TAN ( P’S I 1 ) s V2T=V2A/TAN(PSI2) 

00460 V t =SQRT ( VI A*V I A+V 1 T*V IT ) ? V2-S0RT ( V2 A*V2 A+V2T*V2T ) 

00470 W 1 A=V 1 A ? W2 A=V 2 A ;WIT*V1T-UW tW2T=V2T-UK 
00430 W 1 ~5QRTC W I A*W I A+W 1 T*W 1 T ) SW2=SQRT( Vi 2 A* W2A+K2T+N2T ) 

00490 AZI *SORT( GAMMA+GR AV*RGC*CC)NV*TZl ) .• : 

0500 AZ2“S0RT( GAMMA*GRAV*RGC*C0!W*TZ2 ) 

0051 0 TEMP 1 =TZ f *(l . - (GAMMA- 1 . )* ( V 1 A* VI A+V 1 T*Vi T) / (2 .*AZ1*AZt I) 
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00520 TEMP2=TZ2*( I . -(GAMMA- 1 . >*<V2A*V2A+V2T*V2T>/<2.*AZ2*AZ2>> 

530 530 Al =SQRT ( GAMMA*GR AV*RGC*CONV*TEMP 1 > 
u0540 A2=SQRT ( GAMMA*GRAV*RGG*C0NV*TEMP2 ) 

00550 ARA=A1/A2 nKRA-l ./ARA 
00560 IF( INPUT .EQ. 2) GO TO 580 

00570 Ml A=itl A/A I lM2A=W2A/A2 |MIT=W1T/Al IM2T-W2T/A2 

005 80 580 M I «SQRT CM I A*M1 A+M 1T*M 1 T ) lM2=SQRTCM2A*M2A+M2T*M2T) 

00590 QTC 1 = AT AN2 ( SORT ( 1 ,-M2A*M2A) , -M2A) 

00 600 QTC2 =AT AN2 ( -SORT < 1 . -M2 A*M2A > , -M2 A > 

0061 0 GTC1 *SIN(GTC1 )/(KRA*( 1 .+M2A*C0S(QTC1 )+M2T*SIN(QTCI ) ) ) 

00620 GTC2-SIN( GTC2 )/( KRA*< I .+M2A*C0S(0TC2)+M2T*SIN(QTC2 ) ) > 

00630 GTC=GTC 1 *QTC*QTCI IM*0 
00640 640 M-M+l |QQTC=I 80.*QTC/PI 
00650 ATC.-C 1 .-GTC*tfl T ) * ( 1 . -GTC*MI T)-GTC*GTC 
00 660 BTC*-2 . *GTC*M I A*< 1 , -GTC*M 1 T ) 

00670 CTC=-GTC*GTC*( l .-Ml A*MI A) 

00680 X=BTC*BTC-4.*ATC*CTC 
00690 [F(X .LT. 0.) GO TO 820 

00700 DDTO(-BTC-SIGN< 1 .,BTC)*SQRT(BTC*BTC-4.*ATC*CTC) )/2. 

007 1 0 QC=ATAN2( DDTC, ATC ) t QQC-I 80 . *QC/P I * PR I NT 1 460 
00720 IF<N .EQ. 2) GOTO 740 

00730 PRINT* "NO ZZLE ATTENUATION* 1 ;PRINT 1460 |G0 TO 750 
00740 740 PR INTi “BUCKET AND STAGE ATTENUATION* 1 *PRINT 1460 
00750 750 PR I NT *» UP STREAM CUT-OFF OCCURS AT INCIDENCE ANGLES OF" 
00760 CUT0FF=90.+180.*ARSIN<MJA)/Pt. ' I PRINT 1430, CUTOFF 
00770 CUTOFF=-C UTOFF I PRINT 1430, CUTOFF » PR I NT 1460 
0730 PR INTA“ DOWNSTREAM CUT-OFF OCCURS AT A TRANSMISSION ANGLE OF 
^0790 PRINT 1430, GOTO ? PRINT 14.60 
00800 PR I NT*" THE CORRESPONDING INCIDENCE ANGLE IS" 

00810 PRINT 1430, GQC t PR I NT 1460 

00820 820 IF(M .EQ, 2) GO TO 840 

00830 GTC*GTC2 t QTC*QTC2 t GO TO 640 

00840 840 IF(M1 .LT. !. .ANQ, M2 .LT. 1 .) GO TO 860 

00850 PRINT* "SUPERSONIC RELATIVE FLOW" |G0 TO 1400 

00860 860 IF(N .EQ. 1) GO TO 890 

00870 PR I NT* "THETA THETA T T V R " 

00880 GO TO 900 

00890 890 PRINT*" THETA I THETA T T V R ATTEN" 

00900 900 RHORA«PS 1 *TEMP2/( PS2*JEMP l ) 

00910 ALPHA=ATAN2( Ml T,MI A) I BETA=ATAN2 ( M2T , M2A ) 

00920 GO TO (930, 1056)* N 

00930 930 QQ( J>«90.+180.*ARSIN(M1A)/PI fQQM IN=-OQ( 1 ) |S=IO, 

00940 QQH0LD S «10,*AINT((QQ( 1 > + 10. )/IO. ) 

00950 IFtOQC .LT. 0.) GOTO 980 

00960 QQl't > «-90 .-180 ,-*ARSIN (Ml A ) /PI ?QQMAX=-OQ< I) IS=-1 0. 

00970 OQHOLD-10,*AINT((QQ( 1 1-10. )/10.) 

00980 980 DO 1040 J*2, 100,1 
00990 Q0( J )s*QQHOLD-S |QQHOLD=QQ< J) 

01000 IF(ABS(QQ(J)-QQC) .LT. ABSM.1*S>> S=Q.1*S 
01010 IF(ABS(S) .LT. 0.0001 ) GO TO 1050 
01020 IF( QQC .LT. 0.) GO TO 1040 
01030 IF(OQ(J) .GT. QQMAX ) GO TO 1050 
71 1 040 1040 IF(QQ( J) .LT. QQMIN) GO TO 1050 
01050 1050 DO 1330 K= 1,1 00,1 


ORIGIN^ 
OF 





01060 IF C M .EQ, I) QQSA VH ( K ) =0Q C K ) 

" 070 IF(GOCK) ,LT. QQC .AND. N , HQ. 2) GO TO 1330 
u 1 030 I?(K .HQ. (J+D) GO TO 1340 
01090 Q=Pl*nQ(K)/180. 

01100 IFC AB5( QQCK ) ) .GT. 0.001) GO TO 1110 ;KX0VKY=1 00000. <G() TO 1120 
01110 1110 KXOVKY=( ( I ,+W A*C0S(Q)+M1 T*S I N £ G ) )/(KRA*5IM(Q ) )-M2T)./M2 A 
01120 1120 QR-ATAN C < 1 .—Ml A*M 1 A) *3 IN <Q )/( ( 1 . +M 1 A*M 1 A ) *COS ( Q ) +2 . *K ! A ) > 
01130 G=KRA*STil(0')/ ( I . +W 1 A*COS < 0) +M I T*S I N ( 0 ) ) 

Oil 40 A=( l ,-G*M2T)*( l . -G*M2T )-G*G ; B=-2 . *G*M2A* ( 3 .-G*M2T> 

01 1.50 C=-G*G*( 1 .-M2A*M2A> rDD=C-3-SIGN( 1 . ,B)*SQRT(B*B-4.*A*C ) )/2. 

01160 QT=ATAM2(DD f A) ?0QT=1 G0.*0T/PI 

01170 T1 “ARA*(M2A+C0S(0T) ) iV1=ARA *R1 =C05(QR )-H1 A 5 Cl -?/. 1 A+COSCQ) 

0 1 1 30 T2=RH0RA*( I ,+tf2*C0S (BETA-QT ) ) ? V2=RHORA*CM2A-KXOViCY*jM2T) 

01190 R2=J.U*C0S<ALPHA+0R)-1 . ; C2-1 .+M1 *COS( ALPHA-0 ) 

01200 T3=S i N ( GT -BET A ) l V3=-C SIN C BETA >+i<XOVKY*COS( BETA > ) ?R3=0.;C3=0. 
01210 D= T1 *V2*R3+V| *R2*T3+R UT2*V3-R 1 *V2*T3-TI *R2*V3-V1r*T2*R-3 
01 220 T=(CI *V2*R3+V 1 *R2*C3+R1 *C2*V3-R1 *V2*C3-Cl *R2*V3-V1 *C2*R3) /D 
01230 V=(TI *C2*R3+CI *R2*T3+R1 *T2*C3-R I *C2*T3-T1 *R2*C3-C1 *T2*R3)/D 
01240 S ~<Tt*V2*C3+Vl *C2*T3+C 1 *T2*V3-C I *V2*T3-T 1 *C2*V3-V I *T2*C3 ) /D 
01250 VI A=M t A* A 1 rVlT-M1T*At ;V2A=M2A*A2 ;V2T=M2T*A2 
01260 Z1«(AI +V1 A*C()S ( Q ) +V 1 T*S INC.0 ) j* C A t *COS ( Q ) +V t A ) 

01270 Z2=( A2+V2 A*C05 ( OT j +V2T*S I N ( QT) )*( A2*C()S‘<QT)+V2A) 

0 1 230 AT TEN--4 . 3 42945*AL0G ( T*T*R HOR A*AR A*AR A*AR A* ABS C Z2/Z 1 ) ) 

01290 ATTEN2 CIO =ATTEN2< K) +ATTHK 
01300 GO TO (1310, 1320) ,N 

01310 1310 PRINT 1450, QQC K') ,QQT,T,V,H, ATTEM ?G0 TO 1330 

'1320 1320 PRINT 1440, QQCK) ,QQT,T,V ,R , ATTEN, ATTEN2 CK) ,QQSAVE( K) 

.1330 1330 OQ C K ) = n 0T 

01340 1340 GO TO (1350, !400),N 

01350 3350 IFC INPUT .HO. 2) Go TO 1380 

01360 VI A-V2A |V2A«V3A ?PSI1«PSI2 ;PSI2=PSI3 

01370 P3i~PS2 ?P32-P53 ?TZ1=TZ2 iTZ2=TZ 3 i UN-UN23 ;G0 TO 440 

01330 1330 M1A-M3A ;M2A=M4A ; M1T-M3T «M2T*M4T 

01 390 TEMP 1 -TEMP 3 fT£MP2=TEMP4 ,*PS1=P53 t *pS2=PS4 ,*G0 TO 440 

01400 1400 PRINT 3460 

01410 PR I-ilT i "MOR E INPUT? ( YES«1 , N0=2 ) « ? READ * MORE 
01420 GO TO Cl 50, 1470), MORE 
01430 1 430 FORM AT ( F 1 0 . 3 ) 

01440 1440 FORMAT C2F9. 3, 5F7.3, F9.3) 

01450 1 450 FORMAT ( 2F9 . 3 , 4F7.3) 

01460 1 460 FORMAT ( /) 

01465 111 FORMAT C 1 X , 2 ( F 3 . 3 , IX) ) 

01470 

01430 1470 CONTINUE 
Ol 490 END 
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